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Abstract 

Let TV be a positive integer, c be a positive constant and (U n )n>i be a sequence of 
independent identically distributed pseudo-random variables. We assume that the 
Un's take their values in the discrete set {— TV, —TV + 1, . . . , TV — 1, TV} and that their 
common pseudo-distribution is characterized by the (positive or negative) real num- 
bers 

V{Un = k} = 8m + (-l)* _1 c 

for any k 6 {— TV, —TV + 1, . . . , TV — 1, TV}. Let us finally introduce (S n ) n >o the associ- 
ated pseudo-random walk defined on Z by So = and S n = 53™ =1 Uj for n > 1. 
In this paper, we exhibit some properties of (S n ) n >o- In particular, we explicitly de- 
termine the pseudo-distribution of the first overshooting time of a given threshold for 
{S n )n>o as well as that of the first exit time from a bounded interval. 
Next, with an appropriate normalization, we pass from the pseudo-random walk to 
the pseudo-Brownian motion driven by the high-order heat-type equation d/dt = 
(— l) JV_1 c d 2N /dx 2N . We retrieve the corresponding pseudo-distribution of the first 
overshooting time of a threshold for the pseudo-Brownian motion (Lachal, A.: First 
hitting time and place, monopoles and multipoles for pseudo-processes driven by the 
equation §- t = ±^r. Electron. J. Probab. 12 (2007), 300-353 [MR2299920]). In the 
same way, we get the pseudo-distribution of the first exit time from a bounded inter- 
val for the pseudo-Brownian motion which is a new result for this pseudo-process. 

Keywords: pseudo-random walk ; pseudo-Brownian motion ; first overshooting time ; first exit 
time ; generating function. 
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1 Introduction 

Throughout the paper, we denote by Z the set of integers, by N that of non-negative 
integers and by N* that of positive integers: % = {..., —1, 0, 1, . . . }, N = {0, 1, 2, . . . }, 
N* ={1,2,... }. More generally, for any set of numbers E, we set E* = E\{0}. 

Let TV be a positive integer, c be a positive constant and set k n = (— l)^ -1 . Let 
(LT n )„ S N* be a sequence of independent identically distributed pseudo-random variables 
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taking their values in the set of integers {— JV, — JV + 1, . . . , — 1, 0, 1, . . . , JV — 1, JV}. By 
pseudo-random variable, we mean a measurable function defined on a space endowed 
with a signed measure with a total mass equaling the unity. We assume that the com- 
mon pseudo-distribution of the U n 's is characterized by the (positive or negative) real 
pseudo-probabilities p k = ¥{U n = k} for any k e {-JV, — JV + 1, . . . , JV — 1, JV}. The pa- 
rameters p k sum to the unity: J2k=-NPk = 1- 

Now, let us introduce (S n ) n eN the associated pseudo-random walk defined on Z by 
So = and S n = S?=i Uj f° r n <E W . The infinitesimal generator associated with 
(S n )neN is defined, for any function / defined on Z, as 

fc=-JV 

Here we consider the pseudo-random walk which admits the discrete JV-iterated Lapla- 
cian as a generator infinitesimal. More precisely, by introducing the so-called discrete 
Laplacian A defined, for any function / defined on Z, by 

A/(j) = f(j + 1) - 2f(j) + f(j - 1), j e z, 

the discrete TV-iterated Laplacian is the operator = A o • • ■ o A given by 

N times 

We then choose the pfe's such that Q s = k n cA n which yields, by identification, for any 

k e {-JV, — JV + 1, . . . , —1, 1, . . . , JV — 1, JV}, 

, , k i / 2JV \ /2J^\ 

»=p-* = (-i) "" lc ( fc + iV j and ^ = 1 - C (jvj- (L1) 

When JV = 1, (S n ) n< =^ is the pseudo-random walk to the closest neighbors with a 
possible stay at its current location; it is characterized by the numbers pi = p_\ = c 
and pq = 1 — 2c. Moreover, if < c < 1/2, then p > 0; in this case, we are dealing with 
an ordinary random walk (with positive probabilities). 

Actually, with the additional assumption that p_ fc = p k for any k e {— JV, — JV + 1, . . . , 
JV — 1, JV} (i.e., the U n 's are symmetric, or the pseudo-random walk has no drift), the 
Pk's are the unique numbers such that 

G s f(j) = k n c f^ 2N \j) + terms with higher order derivatives (1.2) 

where f is an analytical extension of f and /( 2Ar ) stands for the (2JV) th derivative of /: 

Our motivation for studying the pseudo-random walk associated with the parameters 
defined by l |l.ip is that it is the discrete counterpart of the pseudo-Brownian motion as 
the classical random walk is for Brownian motion. Let us recall that pseudo-Brownian 
motion is the pseudo-Markov process (X t )t>o with independent and stationary incre- 
ments, associated with the signed heat-type kernel p(t; x) which is the elementary so- 
lution of the high-order heat-type equation d/dt = n N cd 2N /dx 2N . The kernel p(t;x) is 
characterized by its Fourier transform: 



+oo 

E(e luXt ) = I e iux p(t;x) dx = C - Ctu 
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The corresponding infinitesimal generator is given, for any C 2N -function /, by 

G x f(x)= lim hE[f(X h + x)}- f(x)} = n N c^ 2N \x). (1.3) 

h— >0+ tl 

The reader can find an extensive literature on pseudo-Brownian motion. For instance, 
let us guote the works of Beghin, Cammarota, Hochberg, Krylov, Lachal, Nakajima, 
Nikitin, Nishioka, Orsingher, Ragozina, Sato ([1] to [19]) and the references therein. 

We observe that \\-2) and < | 1 . 3 p are closely related to the continuous TV-iterated 
Laplacian d 2N /dx 2N . For N = 2, the operator A 2 is the two-Laplacian related to the 
famous biharmonic functions: in the discrete case, 

A 2 /(j) - fU + 2) - 4/(j + 1) + 6/(j) - 4/(j - 1) + f(j - 2), j G Z, 

and in the continuous case, 

A 2 /W = 0W, xeR. 

In the discrete case, it has been considered by Sato H20H and Vanderbei 02111 . 

The link between the pseudo-random walk and pseudo-Brownian motion is the fol- 
lowing one: when normalizing the pseudo-random walk (S n ) n eN on a grid with small 
spatial step s 2N and temporal step e (i.e., we construct the pseudo-process (eS^/^wj ) {>0 
where [ ■ J denotes the usual floor function), the limiting pseudo-process as e — > + is 
exactly the pseudo-Brownian motion. 

Now, we consider the first overshooting time of a fixed single threshold a < or 
b > (a, b being integers) for (S n )n&i- 

a~ = min{n G N* : S n < a}, a£ = min{n G N* : S n > b} 

as well as the first exit time from a bounded interval (a, b): 

a ab = mm{n G N* : S n < a or S n > b} = min{n G N* : S n <£ (a, b)} 

with the usual convention that min0 = +oo. Since, when at < +oo, S+_, < b — 1 and 

h 

S + > b, the overshoot at time at which is S + — b can take the values 0, 1, 2, . . . , N — 1, 

b b 

that is S a + G {b,b+l,b + 2, . . . ,6+JV-l}. Similarly, whencr; < +oo, S a - G {a-N+l,a- 
N + 2, ...,a}, and when a ab < +oo, S„ ab G {a, a— 1, . . . , a-N+l}{J{b, 6+1, ... ,6+iV-l}. 
We put S+ = S- = ^- and S ab = S aab . 

In the same way, we introduce the first overshooting times of the thresholds a < 
and b > (a, b being now real numbers) for (X f ) t > : 

t- = inf{t > : X ( < a}, t+ = inf{t > : X t > b} 

as well as the first exit time from a bounded interval (a, b): 

r ab = ini{t > : X t < a or X t > b} = mi{t > : X t £ (a, b)} 

with the similar convention that inf = +oo, and we set, when the corresponding time 
is finite, 

X a = X T -, X^ = X T +, X ab = X Tah . 

In this paper we provide a representation for the generating function of the joint 
distributions of the couples (a~,S~), (a^S^) and (a ab . S ab )- In particular, we derive 
simple expressions for the marginal distributions of S~, and S ab . We also obtain 
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explicit expressions for the famous "ruin pseudo-probabilities" P{a~ < cr^} and P{c^ < 
o~a}- The main tool employed in this paper is the use of generating functions. 

Taking the limit as e goes to zero, we retrieve the joint distributions of the couples 
(t~,X~) and (t^~,X^~) obtained in H101 1111 . Therein, we used Spitzer's identity for 
deriving these distributions. Moreover, we obtain the joint distribution of the couple 
(T a b,X a b) which is a new and an important result for the study of pseudo-Brownian 
motion. In particular, we deduce the "ruin pseudo-probabilities" P{r~ < t^} and 
P{tj~ < r~}; the results have been announced without any proof in a survey on pseudo- 
Brownian motion, 111 311 . after a conference held in Madrid (IWAP 2010). 

In HI OH , 01611 and H17B , the authors observed a curious fact concerning the pseudo- 
distributions of X~ and X^ : they are linear combinations of the Dirac distribution and 
its successive derivatives (in the sense of Schwarz distributions). For instance, 

P{X+ G dz}/dz =J2^T S b\ z )- d- 4 ) 

3=0 3 ' 

The quantity 6^ is to be understood as the functional acting on test functions <j> ac- 



cording as (5^',<f>) = (— iy<f>W(b). The appearance of the S^'s in 11.41, which is 



quite surprising for probabilists, can be better understood thanks to the discrete ap- 

(i) 

proach. Indeed, the dj s come from the location at the overshooting time of b for 
the normalized pseudo-random walk: the location takes place in the "cluster" of points 
6,6 + e,6 + 2e, ...,&+ (iV- l)e. 

In order to facilitate the reading of the paper, we have divided it into three parts: 

Part I — Some properties of the pseudo-random walk 

Part II — First overshooting time of a single threshold 



Part III — First exit time from a bounded interval 
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Part I — Some properties of the pseudo-random walk 



2 Pseudo-distribution of Ui and S n 

We consider the pseudo-random walk (S n ) ne n related to a family of real parameters 
{Pk, k e {-N, N}} satisfying p k = p_ k for any k e {1, . . . , N} and J2k=-N Vk = 1. 
Let us recall that the infinitesimal generator associated with (S n ) k > is defined by 



N 



S S M = E p*/0' + fc ) - /C?) = bo - i)/C?) + 5>[/tf + k) + f(j - k)}. 



k=-N 



k = l 



In this section, we look for the values of p k , fee {— N, . . . , N}, for which the infinitesimal 
generator Q is of the form ( |1.2[ ). Next, we provide several properties for the correspond- 
ing pseudo-random walk. 

Suppose that / can be extended into an analytical function /. In this case, we can 
expand 

f(j + k) + f(j-k) = 2j2^f m U)- 



Therefore, 



iV oo 



QM = (po - i)/(i) + 2 E (2^! /(2 ' )(i) 



fc=l ^=0 



JV 



A' 



fe=i 



= 1 \/c=l 



/ (2 %-) 

(2*)! ■ 



Since p + 2J2 k =iPk = we see tri at trie expression |1.2 1 of (/ holds if and only if the 
p k 's satisfy the eguations 

N N 

E fc2 V = for 1 < 1 < N ~ 1 and ^lk 2N p k = -k n c{2N)\. (2.1) 



k=l 



k=l 



Proposition 2.1. The numbers p k , k e {1, . . . , TV}, satisfying \2. l\ are given by 

Pk = (-l) k ~ 1 c[ 

In particular, p N = k n c. 



2N 
k + N 



Proof. First, we recall that the solution of a Vandermonde system of the form J2k=i a k x k 
oil, 1 < t < N, is given by 



Xk 



di, . . . , a N 
ai, . . . ,ajv 



Mai, . . .,a N ) 



1< k < N 



with 



A(ai, . . . , a N ) 



ax . 


■ a N 






aj . 


■ a 2 N 




n N 


af . 


a N 
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and, for any k e {!,... , N}, 



01, . . . , ajv 
ari, . . . ,ajv 



rii 



of 



a Li "2 a 2 k+1 



a fe _ 1 a,v a fe+1 



2 

JV 



(7 



In the notation of and that of forthcoming determinants, we adopt the convention 
that when the index of certain entries in the determinant lies out of the range of k, 
the corresponding column is discarded. That is, for k = 1 and k = N, the respective 
determinants write 



Ol, 
Oil, 



, Ojv 



«1 02 



N 



ajv 



, A 



AT 



Ol, . • . , ffljy 



ai, . 



, aw 



ai 



of 



ajv-i «i 



AT 



It is well-known that, for any k G {1, . . . , AT}, 

A(ai, . . . , ojv) = Y[ a i II ( a ™ - a t) 

l<j<N l<£<m<N 

=(-i) k+N n a i n n 



Kj<N l<j<? 



1<^<TTI<7V 



In the particular case where ag = for 1 < ^ < JV — 1, we have, for any fc e {1, 
that 



,JV}, 





Ol 


dk-1 


a N 


A fe f ai '---' a ^) = (-l) fc+ ^a Ar 
\ai,...,a N J 


a? 


a fc+l 


■ a% 




of" 1 


JV-1 7V-1 
"fc-1 "fe+1 


N-l 
a N 




n °j 

l<j<N 


Yl [a m - a^. 

l<e<m<N 




Therefore, the solution simply writes 








UN 


l< k< N. 





Now, we see that system ( |2.l| l is a Vandermonde system with the choices o^ = k 2 , 
Xk = Pk and ae — for 1 < ^ < JV — 1, — n N c(2N)\/2. With these settings at hands, 
we explicitly have 

k 2 n ( fc -^) n ( k +j) 



ak JJ_ (at- - a,-) = k 2 Y\_ ( k 
1 



2 J 2 ) 



i<i<jv 



l<j<JV 



l<i<jv 



5V* 



(-l)»-*(N + k)\(N-k)\ 



and the result of Proposition |2 . 1 1 ensues . □ 
Finally, the value of p is obtained as follows: by using the fact that J2k=o(~ Ck) = 



0, 



Po 



l- E ^ = l + c E (-!) A 



— N <k< N 
k^O 



1 - c 



/2JV 



- N<k<N 
fc^O 

2iV 

-i)^E(-i) fc 

fc=0 



2JV 



2iV 
k + N 

= 1 - c 



/2JV 
I JV 
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We find it interesting to compute the cumulative sums of the p/s: for k e {—N, . . . , N}, 



E Pi= E 

j=-N j=-N 



5 j0 + 



2N 
3 + N 



k+N 



= l {fe > 0} + (-l) Ar - 1 c^(-iy 



3=0 



'2N 



The last displayed sum is classical and easy to compute by appealing to Pascal's formula 
which leads to a telescopic sum: 



k+N 

E 

3=0 



'2N 
3 



Thus, for k <= {-N, ... ,N}, 

k 



(-1) 



k+N (2N-l 
\k + N )' 



E Pi = Mk>o) + (-1)" - 



j=-JV 



Observe that this sum is nothing but P{£7i < fc}. Next, we compute the total sum of the 
IPj-l's: by using the fact that Y? k H Ck) = 



N 



E 



k=-N 



1 - 



1 - 



(2N 



+ E 



-N<k<N 



2N 

k + N 



(2N 

N 



+ c 



*N 



2N 
N 



= cA N -1 + 2 



1 -c 



As previously, there is an interpretation to this sum: this is the total variation of the 
pseudo-distribution of U\. We can also explicitly determine the generating function of 

Uv for any Ce C*, 



Ht Ul )= E vkC k = i + c 



k=-N 



N 



E (-i)"" 1 



_k=-N 
2N 



2N 

k + N 



= 1 + K„C 



(i-O 



2N 



We sum up below the results we have obtained concerning the pseudo-distribution of 

E/i. 

Proposition 2.2. The pseudo-distribution ofU\ is determined, for k e {-N, . . ., N}, by 

nui = k} = 5 kQ + (-i) k - i c( k m N ] ) 1 

or, equivalently, by 

Wl <k} = l {fe >0} + (- 1 ) fe " lc ( 2 ^ 1 )- 
The total variation of the pseudo-distribution ofU\ is given by 

ri+c[4»-2(»)] if0< c < 

l|Pt " ll = W-i ifcai/ff). 

The generating function ofU\ is given, for any £ e C*, by 

(i-0 2N 



E(C 



1 + K M C- 



(2.2) 
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In particular, the Fourier transform of U± admits the following expression: for any 9 £ 

[0,2*]. by 

E(c iec/l ) =1-04^8^(0/2). 
In the sequel, we shall use the total variation of U\ as an upper bound which we call 

'l + c 



Mi = 



c4 



(2.3) 



4^-2( 2 ^)j ifO<c<l/( 2 ^), 
. N -l ifc>l/( 2 ^). 

Set f(0) = E(c WUl ) for any 9 e [0,2tt]. We notice that f(0) 6 [1 - c4 w ,1] and, more 
precisely, 

ifO<c< l/2 2N -\ 



sup |f (0)| = max(|l- c4 N |,l) 
ee[o,27r] 



1 



cA N -l ifc>l/2 



2N-1 



Let us denote this bound by M^: 



1 ifO<c<l/2 2W - 1 / 
cA N ~l ifc>l/2 2N -\ 



In view of (2.3 i and l2.4l, since 



( 2 #) 



< 2 



2N-1 



we see that Mi > M m > 1. 



Proposition 2.3. The pseudo-distribution of S n is given, for any k e {—Nn, 

by 



e=o 



2N£ 

ej\k + N£ 



(2.4) 



,Nn}, 



(2.5) 



Actually, the foregoing sum is taken over the £ such that I > \k\/N. We also have that 



P{S n < fc} = l{ fc >0} + (-1)* j> 



t.= l 



(2N£ - 1 
£j\k + N£ 



(2.6) 



Proof. By the independence of the U/s which have the same pseudo-probability distri- 
bution, we plainly have that 

E(e ie5 ") =f(0) n = [l-c4 Ar sin 27V (0/2)]". 

Hence, by inverse Fourier transform, we extract that 

f27T 



V{Sn = k} 



1 

2^ 



f(0)™ e - Ifctf d<9 



n 



E(-4 

1=0 



N cY 



n\ 1 
27 



2tt 



sin 2iVf (0/2)e- lfc V0. 



(2.7) 
(2.8) 



By writing sin(0/2) = (c 

,.2tx 



W/2 



A0/2 



) / (2i), we get for the integral lying in ( 2.8 1 that 



1 

2tt Jo 



sin 2iV£ (0/2)e- 



-ikO 



d9 



1 f 2 * fe i9 / 2 -c- ie 



2tt Jo 

t 
(2 



2i 



( -1)™ rvv-iv 



2NI 

e 

2N£ 
m 



-ike 



d0 



j{m-k-Nl)6 d0 



(-1) 



NC 



(-1) 



j{m-k-Nl)6 



m=0 

k / 2N£ 



4 Nt \k + N£ 



(2.9) 
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By plugging ( 2.9 1 into 1 2.8 1, we derive \2.5\ . Next, we write, for k e {— Nn, . . . , Nn}, 
that 



j = -Nn 



kA(Ni) 

E 

j=(-Nn)V(.-Ni) 



2N£ 
j + Nt 



e=o 



k/\{N£) 



2N£ 
j + Nl 



(2.10) 



If k < 0, then the term in sum ( 2.10 1 corresponding to £ = vanishes and 



1=1 v 7 j=-Nl v 7 

The second sum in the foregoing equality is easy to compute: 

y (-iy( me )i= y k-iyf 2A ^ WW 2 ™ -1 ) 



j=-Nl 



= (-l) k ( 2Nl - 1 \ 

K ' \k + N£J 

If k > 0, then the term in sum ( |2.10| l corresponding to £ — is 1 and 



(2.11) 



p{s*<*} = i+£(- C )T; E (- 1 



£=1 



kA(N£) 



j=-Nt 



2N£ 
j + N£ 



By using the convention that (";) = if /3 > a, we see that the second sum above also 
coincides with ( |2.11| l. Formula fl2.6| l ensues in both cases. □ 

Proposition 2.4. The upper bound below holds true: for any positive integer n and any 
integer k, 

\¥{S n = k}\ < y/V{S 2n = 0} < M£. (2.12) 
Assume thatO < c < 1/2 2JV_1 . The asymptotics below holds true-, for any 5 e (0, 1/(2N)), 



F{S n = 0} = O (l/n 5 ) . 
Proof. Let us introduce the usual norms of any suitable function . 



(2.13) 



2tt 



271 



\m\d9, n\\ 2 = 



2n 



2 IT 



\m\ 2 de, u\L= sup 

0G[O,2tt] 



and recall the elementary inequalities 



< UL < 



It is clear from l|2.7[l that, for any integer k, 



\nsn = k}\<~jj\m\ n w=\n 



< n = y 

< iitii: = m, 



i 
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This proves 1 2.12 1. Next, by 1 2.7 1, since f (2tt — 6) = f(6), we have, for any e e (0, it), that 

p{s n = o} = ±- r de = i r f^) n 

f(0) n d0 + / f(0) n d0 



The assumption < c < 1/2 2JV_1 entails that |f(0)| < 1 for any G (0, tt). We see that 
|f < 1 on [Q,e], and |f(0)| < |f(e)| on [e,ir] for any e G (0,tt). Hence, 

|P{«9 n = 0}|< £ + |f( £ )|™. 

Now, choose e = 1/n 5 for a positive 5. We have that 

ln(|f( £ )D = nln(l - c 4 w sin 2A '(^)) ~ -c, 



,1-2AT(5 



which clearly entails, for large enough n, that |f(e)|" < exp (— | n 1 2Ar<5 ). Thus, if £ e 
(0, 1/{2N)), |f(e)| n = o(e) which proves ( [2~13| l. 

Ifc= l/2 2Ar -\ f(0) = l-2sin 2JV (0/2). In this case, the same holds true upon splitting 
the integral ^ into ^ + J^ £ + J*_ e . □ 

Remark 2.5. A better estimate for |P{5„ = fc}| can be obtained in the same way: 

P{S n =0} if n is even, 



Vfc G Z, |P{5„ = fc}| < 



P{5 , n _i = 0} if n is odd. 



Nevertheless, we shall not use it. We also have the following ineguality for the total 
variation of S n : 



Ni. 



E ip{5„ = fc}|<np c/1 ir=Mf. 



k=-Nr. 



Proposition 2.6. For any bounded function F defined on Z r 

\E[F(Si,...,S n )]\ < ||F|LM X ". 



(2.14) 



Proof. Recall that we set pk = V{U\ = k} for any k G {—N, . . . , N}. We extend these 
settings by putting p k = for k G Z\{— iV, . . . , iV}. We have that 



\E[F(S 1 ,...,S n )}\= E ^(^i,...,fc„)P{5 1 = fc 1 ,...,5„ = fc n } 
(fei,...,fc Tl )ez ri 

< II^IL E bfclPfe2-fcl •••Pfen-fcn-ll- 

(fci,...,fc„)GZ' 1 

The foregoing sum can be easily evaluated as follows: 

E \PkiPk 2 -ki ■ ■ ■Pk n -k n - 1 \ 
(fci,...,fc„)GZ™ 



tiGZ \fe 2 ez V \fc„ez 

(Em) =^r 



which proves l |2.14[ l. 



□ 
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3 Generating function of S n 

Let us introduce the generating functions, denned for complex numbers z, by 

G k (z) = p i s n = k}z n = F i S ™ = fc ) z " for k e Z ' 



nGN 

n>|fc|/W 



g(c,z)= ns« = k}c k z n = ns n = k}c k z r ' 



\k\<N n 



We first study the problem of convergence of the foregoing series. We start from 

£ \ns n = k}c k z n \< \c\ k Wooz\ n . 



\k\<Nn 



If C ^ and |C| 7^ 1, then 

ici i-i/ici 



1 / oo oo \ 

^ Ei^r-iciEI M -<T 

\n=0 n=0 / 

- / oo ^ oo \ oo 

+ rrW E - ]Jr E I m»*/cT - E i 

' ISI \n=0 ISI n=0 / n=0 



iMoo^r 

n— C 

If we choose z,C such that IM^I < 1, \M OC) z( N \ < 1 and \M 00 z/C, N \ < 1 (which is 
equivalent to | z\ < ^[min(|C|, l/lCl)]^- or j/\M^z\ < \(\ < l/^f\M^z~\), then the double 
sum defining the function G((, z) is absolutely summable. If \C\ = 1, then 

oo 

ICl' c |Moo^r-^(2iVn+l)|M 00 ^r. 

fceJ,neft n=0 

fc|<JVre 

If we choose z such that |Mooz| < 1, then the same conclusion holds. 
Now, we have that 

G (c^) = E (e p ^ = *K*V = E E (c s "K - E = r-^77ZM 

and, thanks to l |2.2[ l, we can state the following result. 

Proposition 3.1. The double generating function of the ¥{S n = k}, k e Z,n e N, is 

given, for any complex numbers z, ( such that {f\M^z\ < \(\ < l/{f\M^z\, by 



(l-z)C N -k n cz(1-Q 
In particular, for any 9 e [0, 2tt] and z e C such that \z\< l/M, 

1 

1 - z + c4^zsm 2JV (0/2) 

On the other hand, 



G (e ie ' z ) = ; , tM - 2Af ,„,^ - (3.2) 



G (c, *) = E ( E = fc > z " ) c fc = E G ^)c fe - 
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By substituting £ = e l6 in the foregoing eguality, we get the Fourier series of the function 

G(e w ,z)=Y / Gk(z)e ik 



from which we extract the sequence of the coefficients (G k (z)) k ^. Indeed, since P{S n = 
—k} = P{S„ = k}, we have that G k {z) = G- k {z) and 



G k (z) 



2tt 



2tt 



G(e w ,z)e- ik0 d9 = 



2k 



271 



G(e w , z)e lke dO 



2ni 



G{u,z)u k 1 du 



where C is the circle of radius 1 centered at the origin and unclockwise orientated. 
Then, referring to \3.1\ , we obtain, for any zeC satisfying \z\ < that 



G k (z) = 



1 
27ri 



,jk+N-l 
P z (u) 



du 



where P z is the polynomial given by 

P z (u) = (1 - z)u N - k n cz(u - l) 2N . 
We are looking for the roots of P z which lie inside the circle C. For this, we introduce 



the N th roots of k„ : 



'N ■ w 3 



e^N-* for 1 < j < N; 



From now on, in order to simplify the writing of the roots of P z , we make the as- 
sumption that z is a real number lying in (0, l/M^) (and then z e (0, 1)). The roots of P z 
are those of the equations u 2 — 2[1 + 9jw(z)]u + 1, 1 < j < N, where 



w(z) 



T^z 

2Wcz ' 



They can be written as 

Uj(z) = 1 + 8jw(z) — 6jy/w(z) [dj(z) + iejbj(z)], 1 < j < N, 
Vj (z) = 1 + 8jw(z) + 6j tJw(z) [a, (z) + iejbj (z)} , 1 < j < N, 



with 



ej = sgnf sin(^y^7r)^ (with the convention that sgn(0) = 0), 
w(z) 2 - Acos( 2i ^tt)w(z) + 4 + w(z) - 2 cos( 2l ^ 1 ir) 



e 3 


= sgn 


a 3 {z) 


l 




1 

~ V2 



w(z) 2 - Acos( 2j w 1 tt)w(z) + 4 - 10(2) + 2cos( 2j N 1 ir) 



1/2 



1/2 



We notice that a,j(z)bj(z) — | sin (^^7r)|. Because of the last coefficient 1 in the poly- 
nomial u 2 — 2[1 + 0jw(z)]u + 1, it is clear that the roots Uj(z) and Vj{z) are inverse: 

Vj(z) = l/uj{z). 

Let us check that < 1 < for any j e {1, ...,N}. Straightforward 

computations yield that 



\uj(z)\ 2 = Mz) - Bj(z) and \v,(z)\ 2 = A 3 (z) + B 3 (z) 



where 



Aj(z) = w(z) 2 + w(z) yw(z) 2 -4cos(^^7r)w(z) + 4-2cos(^7r) 
Bj(z) = 2^Jw(z) (w(z) — cos(^^7r))a :) -(z) + | sin(^^7r) | bj(z) 
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Since Vj(z) = l/u (z), checking that |i*j(2r)| < 1 < |uj(-z)| is equivalent to checking 
that \uj{z)\ 2 < \vj{z)\ 2 , that is B 3 (z) > 0. If sin(^Tr) ^ 0, we have aj{z)bj{z) ^ 0, 
cij(z) = | sin(^^7r) \/bj(z) and then 

Bj(z) = 2aj(z)y/w(z) w(z) - cos(^^tt) + bj(z) 2 

= a J (z)y / w(z) \Jw(z) 2 - 4cos(^i7r)w(z) + 4 + w(z) > 0. 

If sin(^=^7r) = (which happens only when TV is odd and j = (N + l)/2), a (z) = 
^Jw(z) + 2, bj(z) = and then Bj = 2a j (z) v / w(z)(w(z) + 1) > 0. 

The above discussion ensures that the roots we are looking for (i.e. those lying inside 
C) are Uj(z), 1 < j < N; we discard the Vj(z)'s. 
Remark 3.2. We notice that 



lim dj(z) — \/2sin(^^7r) , lim bj(z) = tj\f2 cos(^-7r 

— s-1- z— >\~ 



Z-S-l 



and then 



lim 9j [ aj (z) + \t 3 b 3 {z)] = V^ifj 



where we set <pj = -ic'Tjr 1 . The ipj, 1 < j < N are the (2N) th roots of k n with positive 
real part: ip 2N = k n and $l(tpj) > 0. As a result, we derive the asymptotic, which will be 
used further, 



i j {z) = l + s j {z) with Sj(z) 



2 VT 



Example 3.3. For N — 2, the roots explicitly write as 

u\{z) = [1 — bi(z)w(z)] — i w(z) — a\(z)y/w{z) 



(3.3) 



u 2 {z)=ui(z), vi(z) = l/ui(z), v 2 (z)=vi(z), 



with 



w{z) = ^J-, ai (z) = ^^w(z) 2 +4 + w(7), b 1 {z) = ^^Jw{z) 2 +A-w{ 



2y/cz ' iV ~ y y/2 
For N = 3, the roots explicitly write as 

1 r 



ui{z) = ^ [2 - w{z) - (ai(z) - 6i(«)v^) V^U) 

- 1 - \w{z)Vl- (ai(«)>/3 + 6i(«)) 
u 2 (z) = 1 - w(z) - a 2 {z)y/w(z), u 3 (z) = u\(z), 



with 



vi(z) = l/ui(z), v 2 (z) = l/u 2 (z), v 3 (z)=vi(z), 



ai{z) = ^j=\J yjw{z) 2 - 2w{z) +i + w{z) - 1, 
h(z) = ^j=\j yjw{z) 2 - 2w{z) + 4 - w{z) + 1, 
a 2 (z) = v / w(z) + 2, w{z) = ■ 
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Now, Gk(z) can be evaluated by residues theorem. Suppose first that k > (then 
k + N — 1 > 0) so that is not a pole in the integral defining Gk(z): 



' P z(u) J <H P' z {Uj{z)) 



3 = 1 

N 

I II ; I 7. I 

u 2 {zf. (3.4) 



1 ^ l-u 3 (z) , 
N{l-z)j^l + u j {z) ' 



The foregoing representation of Gk{z) is valid a priori for any z e (0, l/M^). Actually 
in view of the expressions of w(z) and Uj(z), we can see that (3.4 1 defines an analytical 
function in the interval (0, 1). Since Gk(z) is a power series, by analytical continuation, 
equality ( 3.4 1 holds true for any z e (0, 1). Moreover, by symmetry, we have that Gk{z) = 
G-k{z) for k < 0. We display this result in the theorem below. 

Theorem 3.4. For any k e Z, the generating function of the ¥{S n = k}, n e N, is given, 
for any z € (0, 1), by 

G k {z) = ^44%W |fe| - (3-5) 

N(l -z) ^ 1 + Uj{z) 



Remark 3.5. Another proof of Theorem |3 .4| consists in expanding the rational fraction 
C i-> G(£, z) into partial fractions. We find it interesting to outline the main steps of this 
method. We can write that 

<*<.*> -E^+E B ' iz) 



with 



Uj{z) l-Uj{z) ,. = Vjjz) l-Vjjz) = l/uj(z) l-Uj(z) 

''''' N(l- z) l + Uj{zY 3 "' N(l-z)l + Vj(z) N(l-z) 1 + Uj{z)' 



We next expand the partial fractions ^ and <^_J.( Z ) into power series as follows. 
We have checked that \uj(z)\ < 1 < \v/(z)\ for any j e {1, ... ,N}. Now, if \uj(z)\ < |C| < 
\vj(z)\ for any j <= {!,..., AT}, 



from which (|3.5[l can be easily extracted. 



4 Limiting pseudo-process 

In this section, by pseudo-process it is meant a continuous-time process driven by a 
signed measure. Actually, this object is not properly defined on all continuous times but 
only on dyadic times k/2\ j, k e N. A proper definition consists in seeing it as the limit 
of a step process associated with the observations of the pseudo-process on the dyadic 
times. We refer the reader to [11] and [17] for precise details which are cumbersome 
to reproduce here. 

Below, we give an ad hoc definition for the convergence of a family of pseudo- 
processes ((X%) t >o)e>o towards a pseudo-process (X t ) t > . 

Definition 4.1. Let ((X t £ ) t > ) e >o be a family of pseudo-processes and (X t ) t >o be a 
pseudo-process. We say that 

(Xt)t>o — ► (X t )t>o 
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if and only if 

Vn G N*, Vii, ...,t n > 0, V/ii, ...,/i n G K, E^ 2 ^ 1 Mfc ^) — » E(V^=i wXi * ) . 

This is the weak convergence of the finite-dimensional projections of the family of 
pseudo-processes. 

In this part, we choose for the family ((X £ ) t > ) e >o the continuous-time pseudo- 
processes defined, for any e > 0, by 

X t = £ Slt/e 2N \ ) t > 0, 

where [ • J stands for the usual floor function. The quantity Xf takes its values on 
the discrete set el. Roughly speaking, we normalize the pseudo-random walk on the 
timex space grid e 2N NxeZ. Let (X t ) t >o be the pseudo-Brownian motion. It is character- 
ized by the following property: for any n G N*, any t\, . . . ,t n > such that t\ < ■■■ < t n 
and any fxi, . . . , /z n G M, 

E^ e i£fc=i flkXt k^j = e - c K?=i(Mi+---+A'fc) 2JY (tfc-tfc-i) > (4i) 

We refer to 11110 and 01711 for a proper definition of pseudo-Brownian motion, and to 
references therein for interesting properties of this pseudo-process. 

Theorem 4.2. Suppose thatc < 1/2 2JV_1 . The following convergence holds: 

(XI) t >o — ► (X t ) t >o- 

Proof. • We begin by computing the Laplace-Fourier transform of Xf. By definition of 
Xf, we have that E(c^ x ') = E(e iMeS L*A 2N J \ and then 



J e" At E(e^) dt = [ J e " At dt J E(e^ s ") 

oo 

^(e-^ 2 ")"E(e i ^) 

n=0 

£ (e-^ N )>-f P{S„ = fc} 
A G(e^,e-^ 2N ). 



n=0 

1 - e~ Ae 
A 

l - e~ Ae 



I -c 



-As 2 



By ( ]3.2p , we have that 

V 6 ' 6 J l-c-^ 2Ar +c4^e-^ 2N sin 2JV (//e/2)" 



Actually, equality \A.2\ is valid for A such that e~ A£ < that is, A > (lnM 00 )/e 2N . 

Since c is assumed not to be greater than 1/2 2N ~ 1 , by ( 2.4 1, we have that Afoo = 1 and 
(|4.2p is valid for any A > 0. 



_2JV 



Now, by using the elementary asymptotics sm(pe/2) = [ie/2+o(e) and e Ae - 

E->0+ E-S-0 + 

1 — Ae 2Ar + o(e), we obtain that 
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As a result, for any A > 0, 

r+oo -i /-+00 

lim / e" At E(e* x 'A dt = ^ = / e^^' dt 



from which and l |4.1[ l we deduce that 

lim E(e^ x ') = e^ 2 ™ 4 = E(e^ Xt ). (4.3) 
Notice that the Laplace-Fourier of X t takes the simple form 

r + oo i 

/ e- At E(e i ^)dt= 2N . 
Jo \ + c^ N 

• In the same way we compute the Laplace-Fourier transform of X^ +e2n which will 

be used further. We have E ^e ifJ " x *+ E ' 2 " \ = E(e iA1£S W e 2N j+i) . Then 

J e- At E(e 1 ^+^"J dt = ( J e- M dt\ E(e i ^ 5 "+ 1 ) 



1 - e" 



A 

1 - C Ae X£ 2N ^ - < , J \ x " 



= eA ^- 1 [G(e i ^e-^)-l . 
As for ( |4.3p , we immediately extract the following limit: 

lim_ E ( e i/lX '+« 2JV ) = E (e^ Xt ) . (4.4) 

• We now compute the joint Fourier transform of Xf 2 ) for two times ti,t 2 such 
that t\ <t 2 . Using the elementary fact that [x\ — [y] € { [x — y\ , [a; — y\ + 1}, we observe 
that 

X t 2 ~ X ti = eS Yt2/e 2N \-\_t 1 /e 2N \ € {e5'L(t2-ti)/e 2N J ' e5, L(t2-ti)/e 2JV J+l} = ' ^"ta-ti+e^}' 

Then, we get, for /xi,/X2 € K, that 

By < |4 . 3p and ( |4.4p , we obtain the following limit: 

lim E(e iM2X *2-*i) = lim + E(e iAI2 ^-ti+ e 2 ~) = E ( e i ^ x *2-ti ) 

which yields that 

lim Efe 1 ^ 1 *^ 2 ***)) = lim fife^ 1 ^ 2 ^ ") x lim E(e iAI2A V*i ) 

• Finally, we can easily extend the foregoing limiting result by recurrence as follows: 
for n G W, m, . . . , /i n S M and for any times ti, . . . , t„ such that t\ < ■ ■ ■ < t n , 

lim E^O^+'-'+^O) =E (c i ^ x ' 1 +-+^ x '")). 

The proof of Theorem |4.2] is complete. □ 
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We find it interesting to compute in a similar way the A-potential of the pseudo- 
process (X t )t>o- By definition of Xf, we have, for any ot,/3 € K such that a < ft, 

F{Xf £ [a, (3)} = P{S lt/£ 2 N] e [a/e, 0/e)}. Thus, 

oo / ^(ii+l)^ 2 



/ e - At P{Xfe[a,/3)}di = £ / e" A *dt 



<■ MtJ P{S„ G [a/e,0/e)} 



1 - e 



1-e" 



-Ae 2 



oo 

E (•-»■") 



-Ae^ 



E 

fceZ: 
*/e<fc</3/e 



F i s * = k } 

/e<fc</3/e 

^F{5 n = fc }(e-- 2N )" 



1-e 



-Xe 2 



x/e<k<p/e 



Interchanging the two sums in the above computations is justified by the fact that the 
series X^^Lo^i^™ = k} le~ Xe j is absolutely convergent because of the condition 



c < 1/2 2N ~ 1 . Indeed, by \2.12) , for any A > 0, |P{5„ = k}\ < M™ = 1. 
Put Uj(A, e) = uj ^e~ Ae2N j. This yields that 



N 



Suppose, e.g., that < a < (3. Then, 



3=1 n ' ^ a/e<k<P/s 



a/e<k<p/e k=[a/e~} 



k _ uj( \, e)' a ' E ' —Uj(\,e) 
1 — itj(A, e) 



where [ . ] stands for the usual ceiling function. By using ( 3.3 1, we deduce that 



Uj(X,e) = 1 — v/v X/ce + o(e) 

e->0+ 



(4.5) 



which implies that 
Therefore, 



lim Uj(X,e) 



lim 

e^0+ 



N 



,-At 



3 = 1 

1 N ( c fi 

2 =rE(^/- 



The case a < /3 < is similar to treat. We have obtained the following result. 
Proposition 4.3. The X-potential of the pseudo-process (X t )t>o is given by 



r+oo 

/ e _At (¥{X t € dx}/dx) At = < 
Jo 



1 V w . e -^^ ifx>0 

Vl-l/(2iV) Z^J e " X - U ' 



27VT/^A 1 - 1 /(2A') 



3=1 



1 V^-eK^ 1 ifx<0 



2iV\Y^A 1 - 1 /(2A r > 



3=1 
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Part II — First overshooting time of a single threshold 



5 On the pseudo-distribution of (at, S, 



b > u b 

Let b £ N*. In this section, we explicitly compute the generating function of (a^ , S^). 



Set, for I € {b, b + 1, . . . , b + N - 1}, 

H+ £ (z) = e(^ + 1 {s ^ +<+oo} ) = X>K + = k,S+ = £}z k . 



feGN 

We are able to provide an explicit expression of H^ t (z). Before tackling this problem, 
we need an a priori estimate for P{u^ = k,S£ = £}■ By 1 2. 14 1, we immediately derive 
that \P{a+ = k, 5+ = l}\ = |P{S*i < b, . . . , S k -i < b, S k = l}\ < M*. Hence, the power 
series defining H^ t (z) absolutely converges for \z\ < 1/M\. 

5.1 Joint pseudo-distribution of (a^,S^) 

Theorem 5.1. The pseudo-distribution of(a£, S£) is characterized by the identity, valid 
For any z e (0,1) and any £ e {b, b + 1, . . . , N - 1}, 



E 



where s+ (z) = 1 and for k € {1, . . . , N}, I € {1, . . . , N - 1}, 

S fc,f( z ) = E u h( z )'-- u it( z )i Pk( z )= II i u k(^)-Uj(z)]. 

l<i 1 <--<ip<N l<j<N 
i].,. j^k 

Proof. Pick an integer k > b. If S n = k, then an overshoot of the threshold b occurs 
before time n: < n. This remark and the independence of the increments of the 
pseudo-random walk entail that 

n b+N-1 

F{S n = k}=P{S n = k,a+ < n } = Y, E HS n = k,a+ = j,S+ = £} 

j=0 i=b 

n b+N-1 

= E E n°i=3,st=i}ns n -i = k-i}. (5.2) 

j=o e=b 

Since the series defining Gu{z) and H^ t (z) absolutely converge respectively for z £ 
(0,1) and \z\ < 1/Mi, and since Mi > 1, we can apply the generating function to the 
convolution equality l |5.2[ ). We get, for z e (0, 1/Mi), that 

6+JV-l 

G k (z)= £ Gk-e(z)H+ e (z). 
e=b 



Using expression 13.5 1 of Gu, namely G k (z) = J2 J= i a j( z ) u j( z ) k f° r k >0, where aj(z) 

W{hr) ttstM- we obtain that 

N /b+N-1 Tj+ 

H b,i 

Uj (z) 



J2aj(z)uj{z) k ( ( V - 1 = °> k> b + N-1. (5.3) 



j=i \ e=b 
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Recalling that Vj(z) = l/uj(z) and setting a (z) = oij{z)\Y^^ 1 H^ e (z)v 3 (z) e - lj, 
system (5.3 1 reads J2f=i ®j( z ) u j( z ) k = 0, k>b + N — 1. When limiting the range of k to 



the set 



b + N, b + N + 1, . . . , b + 2N — 1}, this becomes an homogeneous Vandermonde 
system whose solution is trivial: &j(z) = 0, 1 < j < N. Thus, we get the Vandermonde 
system below: 

b+N-1 

£ H+t = 1, l<j<N. 



(5.4) 



System ( |5.4[ l can be explicitly solved. In order to simplify the settings, we shall omit the 
variable z in the sequel of the proof. It is convenient to rewrite ( |5.4[ l as 



b+N-l 



]T ff+«j- 6 = 4 l<j<N. 



l=b 



Cramer's formulae yield 



+ = Vejvi, ...,v N ) 
b ^ V(vi,...,v N ) 



, b<e<b+N -l 



(5.5) 



(5.6) 



where 



V(vi, . . .,v N ) = 



1 vi 

1 t>2 



1 V N 

and, for any I e {b, . . . , b + N - 1}, 



..N-l 



"AT 



n - ^ 

l<i<j<N 



Vl(vi, . . . ,V N ) = 



1 Wl 
1 «2 

1 Wat 



S— 6— 1 



1 w l v l 

i~ b -i u b vt b+1 



l-b+l 



u N u N u N 



..N-l 



..N-l 



This last determinant can be expanded as Ylk=i u k^t{i>i, 
for ke {!,..., N}, 



... v N 

,Vk-l,Vk+l, 



,vn) with, 



Vu{v\,. ■ . , Vk-i, Vk+i, • ■ ■ , v N ) 



1 


Vl 




1 







-1 




1 


Vk-1 ■ 


■ 


1 







-1 


■ 


1 




■ <~ b - 


1 


1 




-1 




1 


Vk+l ■ 


u k+l 


1 





u k+l 


-1 


u k+l 


1 


VN 


J.-h- 
U N 


1 





J.-M 
N 


-1 


N 



In fact, the quantity Vm(vi, . . . , v k -i,v k +i, ■ . ■ , vn) is the coefficient of x l b in the poly- 
nomial 

1 vt 



,N-1 



1 Vk-1 

1 X 

1 v k+1 

1 V N 



N-l 
J k-1 
V N-1 

N-l 
J k+1 



N-l 

J N 
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which is nothing but V(vi, . . . , v^-i, x, ujt+i, . . . , vn), the value of which is 

n ~ v ^ n ^ ~ ^) n ^ ~ ^ 



\<i<j<N 



l<i<k-l k+l<i<N 
Ill<i<j<jv( w i-«i) 



j I (af-wO = V(«i,...,«Ar) JJ 



itk 



l<i<N 
i^k 



X — Vi 



Vk - Vi 



l<i<N 
i^k 



Uk 



UiX — 1 

u k - Ui 



(-1) 



N-l k 



u^ 1 



Pk 



V(vi,...,v N ) Yi {uiX-l). 



l<i<N 

i^k 



Using the elementary expansion Y\i<i< N (mx — 1) = Sfco 1 (~^) A ' 1 is ti xi > we obtain 
by identification that 



Vm{v\, . . .,v k - 1 ,v k +%, ■ • -,v N ) = (-1)* 



b ^ zr s+ t _ b V(vi,...,v N ). 
Pk 



Plugging this expression into 1 5.6 1, we then derive for H^ f {z) representation 1 5.1 1 which 
is valid at least for z e (0, 1/Mx). Finally we observe that (5.1l defines an analytical 
function in (0, 1) and that H^ e (z) is a power series. Thus, by analytical continuation, 
( |5.1| l holds true for any z e (0, 1). □ 

Example 5.2. For N = 2, the settings of Theorem |5 . 1 1 write 

s t,o( z ) = s to( z ) = !> s t,i( z ) = M z ), 4A Z ) = M z ), 
pf(z) = Ux(z) - u 2 (z), p 2 (z) = u 2 (z) - Ui(z), 



where ui{z) and u 2 (z) are given in Example 3.3 and (5.1 1 reads 



\S+=b,< r +<+ 0O } 



3( " , . T 



IL {S+=fc+l,<T+<+00} 



Ul (z) b+1 ~u 2 {z) b+l 

ui(z) - u 2 (z) 
u 1 (z)u 2 (z) b+1 - u 2 {z)u 1 {z) b+1 
u\(z) — u 2 (z) 



Remark 5.3. We have the similar expression related to a a below. The analogous system 
to l |5.4[ l writes as 

a 

e-a-N+l 

The solution is given by 



where H a i (z) = E(^l {s - = ^ CT - <+oo} 



E z a - 1 



{S~ =£,a-<+oo) 



N <T 



-«(*) 



u k {z) 



a+N-1 



where s~ (z) = 1 and, for k € {1, . . . , N}, I € {1, . . . ,N - 1}, 

s k,t( z ) = £ v il (z)---v lf (z) 1 Pk(z)= I [ [Ufc(z) - 



l<i 1 <-.-<i 4 <JV 



l<j<N 



From pseudo-random walk to pseudo-Brownian motion 



21 



The double generating function of (a^,S^) denned by 

b+N-l 



E(z a » ( s » l {(T + <+oo} ) = E E K 61 



admits an interesting representation by means of Lagrange interpolation polynomials 
that we display in the theorem below. 

Theorem 5.4. The double generating function of(er^, S£) is given, for any z e (0, 1/Mi) 
and C G C, by 



E 



(^ + C s >%+ <+oo} ) = X>(z,C)M*)C)* 



(5.7) 



where 



L fe (z,o= n 



k=X 

J 3 < N v k{z)-Vj{zY 



ke {!,... ,N}, 



are the Lagrange interpolation polynomials with respect to the variable £ such that 

L k (z,Vj(z)) = S jk . 

Proof. By ( |5.6[ ) and by omitting the variable z as previously we have that 

b+N-X b+N-l 

e( : 



^V>% b+<+00} ) = E K^ l = E ^'-'^V 



6+JV-l / JV 



E E4 



i=b \ti Vfa, 
N /b+N-l 



E E 

k=X \ £=6 

TV 



it? r c 

V{vi,...,vn) / 



E 

fc=l 



V^l , ■ ■ ■ ; ; Ci Vk+X , ■ ■ ■ , V N ) 

V(vi, . . .,v N ) 



(u k () b 



It is clear that the guantity V(vi, . . . , v k _i, (, v k+ i, . . . , wat)/1 / (wi, . . . , v N ), which explic- 
itly writes as 



1 v x 

1 Wfc-1 

i c 

1 



N-X 



N-X 
J k-X 
rN-X 

>-l 
J fe+1 



AT-1 



1 «1 

1 v N 



N-X 



N-X 



J N 



defines a polynomial of the variable ( of degree N — 1 which vanishes at v\, . . . , v k -x, 
Vk+xt ■ ■ ■ ,vn and eguals 1 atv k . Hence, by putting back the variable z, it coincides with 
the Lagrange polynomial L k (z, Q and formula 1 5.7 1 immediately ensues. □ 

Example 5.5. For N = 2, ( |5.7| ) reads 

= ,f ( A ( Ul (z) b ^- u 2 (z)^) 
ui(z) - u 2 {z) 

+ (u 1 (z)u 2 (z) b+1 -u 2 (z)u 1 (z) b+1 )C]. 
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This is in good agreement with the formulae of Example [5?2| We retrieve a result of [20 1. 



5.2 Pseudo-distribution of 

In order to derive the pseudo-distribution of which is characterized by the num- 
bers Hfr £ (l),£ € {b, 6+1, . . . , b + N — 1}, we solve the system obtained by taking the limit 
in (5^1 as z — >• 1~. 

Lemma 5.6. The following system holds: 

\H+(l) = (-l) k 



b+N-l 

E 

i=k+b 



i-b 
k 



b + k-1 
b-l 



0<k<N-l. 



Proof. By d3.3 1, we have the expansion Vj(z) — l/uj(z) = 1 — £j(z) where £j(z) 



(5.8) 



z-yl- 



<D( 2 tyl - z) for any j e {1, . . .,N}. Putting this into 15.41, we get that 



b+N-l 



£ (l-e j (z)) i -»H+ i (z) = (l-e j (z))- b , 



that is, 



N-l 



(=b 



'b+N-l 



EHH E 



fe=0 



i=b+k 



i-b 



fc=0 



E 

fe=0 



-b 



b + k-l 

b-l 



Set 



b+N-l 

M k (z) = (-l) k 

£=b+k 



i-b 
k 



b + k-1 
b-l 



e=N 



ej{z) k 



b + i-1 
b-l 



(5.9) 



e 3 {zf. 



Then, eguality ( |5.9[ l reads 

w-i 

J2 M k (z)e 3 (z) k = R 3 {z), 1<j<N. 

k=0 

This is a Vandermonde system the solution of which is given by 

Vkl 



M k (z) 



£i(z),...,e N (z) 
R 1 (z),...,R N (z) 

y(ei(z),...,e JV (z)) 



< k < N - 1 



where 









1 ei(«) 


... ei{z) N ~ l 












1 e 2 («) 


... e 2 {z) N - 1 








• -,e N (z)) = 


1 










1 ejv(^) 


... e N (z) N ~ l 








1 


£l(z) ... 




Ri(z) e 1 (z) k 


+1 




,£jv(z) \ 


1 


£2(2) ... 


e 2 (z) k - 1 


R2(z) e 2 (z) k + 1 . 


■ e 2 {z) N - 1 


,R N {z)J 
















1 


e N {z) ... 


SNizf- 1 


R N (z) e N (z) k+1 . 


■ e N {z) N ^ 
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Since, by (3.3 1, Sj(z) ~ constant x 2 {^1 — z for any j e {!,••■ , N}, we have that 



z->l- 



V , (e 1 (z),...,e JV (z))= TT [e m (*) -£*(«)] - constant x (1 - z)^' 1 )/ 4 
and second, 



l<£<m<Ar 



. , /6 + iV-l\ , Njv , 

Rj{z) ~ I j£j(z) ^ constant x y 1 — z 

which implies, for fee {0, . . . , N — 1}, that 

Vk( u l( r Z l , '" ,6 u { f\) = 0\(l-z) l ' 2+ ^ m <N-^ k rn)/{2N)-\ = Q l (l _ 

\Hi(z), . . . ,Kn{z)J z^i- L J L 

Therefore lim z _ i . 1 - Mk{z) — 0. On the other hand, for z € (0,1), referring to the def- 
inition of M k (z), we can see that the quantity H^ e (z) can be expressed as a linear 
combination of the Mk(z)'s plus a constant. Hence, the limit \im z ^ 1 - H^ t {z) exists and, 
by appealing to a Tauberian theorem, it coincides with H^ t {\). This finishes the proof 

of <[5~5T>. □ 

Theorem 5.7. The pseudo-distribution of is characterized by the following pseudo- 
probabilities: for any £ 6 {b, b + 1, . . . , b + N — 1}, 



b fN - 1\ (b + N - 1 



P{5+ =l,a+< +^} = (-l) b+i £{ i _ J) 
Moreover, P{<r^ < +00} = 1. 

Proof. We explicitly solve system ( |5.8| l rewritten as 



(5.10) 



6-1 



as an 



The matrix of the system is ( (f) ) which admits ( (— l) k + e (f) 

/ o<k,e<N-i V yk/ / o<k,e<N-i 

inverse with the convention of settings (. ) = if k > £. The solution of the system is 

given, for I e {0, 1, . . . , N - 1}, by 



N-l /.x 

^Wi) = E(- 1 ) fc+ \^ x (- 1 ) fe 



b+k-i\ . .jb+i-i\^ (b+k-i 

0-1 > ={ - 1} { »-i J£U<-i 



=(-iri' i>+ '" 1 V i,+w " 1N ) =(-i)'— f w - i Y» +w - 1 



This proves 1 5.10 1. Now, by summing the V{S^ = £, <r+ < +00}, b < £ <b + N -1, given 



by ( |5.10[ ), we obtain that 



7 fcb v 



iV-1 



6+iV-lVyv (-1)^ fiV-1 



£=0 
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Writing ^ = f Q x e+b 1 dx, we see that 



JV-l 

E 



(-l) e (N - 1 



l + b V £ 



(l-if-V-'di 



5.7 



x b ~ x dx 

(N-l)\(b-l)\ 
{b + N -1)1 ' 

is finished. 



Hence P{ct^ < +00} = 1. The proof of Theorem 

In the sequel, when considering S£, we shall omit the condition a£ < +00. 
Example 5.8. Let us have a look on the particular values 1, 2, 3, 4 of N. 

• Case N = 1. Evidently, in this case = b and then 

P{S+ = b} = l. 
This is the case of the ordinary random walk! 

• Case N = 2. In this case S£ £ {b, b + 1} and 

H$b =b} = b+l, ¥{S+ = b+l} = -b. 

• Case N = 3. In this case S£ £ {b, b + 1, b + 2} and 

P{S+ = 6} = i(6+l)(6+2), P{5+ = 6+1} = -6(6+2), P{5+ = 6+2} = j 

• Case TV = 4. In this case € {6, 6 + 1, 6 + 2, 6 + 3} and 



□ 



6(6+1) 



P{S^ = b} = l(b+ 1)(6 + 2)(6 + 3), P{5+ = 6 + 1} = -1 6(6 + 2)(6 + 3), 
P{5+ = 6 + 2} = i 6(6 + 1)(6 + 3), P{5+ = 6 + 3} = ~ 6(6 + 1)(6 + 2). 

5.3 Pseudo-moments of 

In the sequel, we use the notation (i) n = i(i — l)(i — 2) • • • (i — n + 1) for any i £Z and 
any n £ N* and (i)o = 1. Of course, (i) n — — n)! and (i) n /n\ = (*) if i > n. We also 
use the conventions 1/i! = for any negative integer i and J2i=i = if i > j. 

In this section, we compute several functionals related to the pseudo-moments of S£. 
More precisely, we provide formulae for E[(5^ — £)1 (Theorem 5.101, E[(S^ — 6) J 
(Corollary[5~Tll, E[(S£) ] and E[(S' fc + )"] (Theorem"" 



5.121. 



Putting the elementary identity 1/(1 + 6) = x e+b 1 da; into the equality 



b + N -1 





b+N-l 



E (-D 



b/N-1 



b + N 
b 



l=b 
N-l 



i\£-b 



f(i) 



we get the following integral representation of E[/ (S^)] . 
Theorem 5.9. For any function f defined on {6, . . . , 6 + N — 1}, 



N - 1 

I 



l(l + b)x l x"' 1 dx. 



(5.11) 
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Theorem 5.10. For any integers n > and ft, the factorial pseudo-moment of (St — 0) 
of order n is given by 

P-^SW + -»)lU> 

r-iv /n\ 

53 (* + 6 - 1)109 - 6 + »- Jb-l)!( ) if (3 > 6 + 1. 



= < 



(6- !)!(/?- 6-1)! ^ ^ '\k 



If n < N — 1, we simply have that 

E [( S b ~ P)J = (-P)n = (-1) W + 1) • • • 09 + n - 1). (5.13) 
Proof. By ( |5.11[ ), we have that 

E[(tf - /9)J = 6( & + N b ' J ) jT 1 (J>1)' 7 ') (/ + * " /9)» ^ +6 - 1 ) dx. (5.14) 

Next, by observing that (£ + b- (3) n x l+b ~ x = x n+l3 - x £^ (x e+b ^ 13 ), we obtain that 

ec- 1 )' r 7 o « + 6 - a- ^ = e r 7 ') ^ (^ +b ^) x 

£=0 V 7 £=0 V V / 



d™ . , / A' - i 

d.r 

d 



E(-^ 

, -((l-x^-V-^)^- 1 . (5.15) 



Applying Leibniz rule to ( |5.15| l, we see that 

1=0 ^ 7 

d n-fc 



dx™ fc 

nA(iV— 1) 

E (-!)* (I) ( N - ~ ^ W 1 - a; ) iV " 1 " 

fe=0 



Therefore, 



1 /JV-1 

Jo \e=o 



N - 1 



+ 6- / 3)„x £+b - 1 I da; 
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nA(N- 


1) 


E 


(-lf(N- 


k=0 




n/\(N- 


1) 


E 


(-l) fc (iV- 


k=0 




'(6- 





(l-x) 



N-l-kk+b-1 



x K+ °- L dx 



nA(JV-l) 



E (- 1 )' 



(fe + 6-1)! fn 



if P < b, 



= < 



(6 + iV-l)! ^ v ' (k + b- 8-n)\\k. 

f-lW-lV " A ^ _1) AA 

fIi^ £ (* + 6- w -6 + „-*-i)ig i^>6 + i. 

(5.16) 



Finally, plugging l |5.16[ l into ( |5.15| l and l |5.14[ l yields l |5.12[ l. 

Assume now that n < N — 1 and j3 < b. If n > 1 — [3, we can write in I 5.12| l that 



(it + 6-1)! d' 3+r 



„k+b-l\ 



(k + b- P -n)\ dx?+ r 



Then, 



nA(JV— 1) 

E (-D 

fc=0V(n+/3-&) 



k (k + b— 1)! fn 



EC- 1 ) 



(fc + 6-/3-n)! Vfc 
(fe + 6-1)! /n 



(fc + 6-/3-n)! Vfc 



E(-d a 



vfc=0 



J/3+n-l 



n — k^.k+b—ft—n') 



k=0 

2(-l)*(n) fc (6-l)^ +n _ fc _i 

(- i rWn(^-D^( /3+ r 1 )=(-l) (6-0)1(0 -1)1 
(-l)"|^|w + l)---(^ + n-l). 



^(l-^V- 1 ) 



/3 + 71- 1 
fc 



)((!-*) 
„(6-l)!(/3 + n-l)! 



(5.17) 



Putting < |5~T7| l into | [5~T2| yields j5TT3) . If n < -/3 (which reguires that (3 < 0), in $5A2\ , 
we write instead that 



(fc + 6-1)! _ 1 

(A: + b-/3-n)! _ (-n-/3)! ./„ 



c fc+h - 1 (l-a;)-"-^d2;. 



Then, 



nA(JV-l) 



E (-D 



fc (fc + 6-1)! Ai 



fc=0V(n+/3-b) 



(fc + 6-/9-n)! \fc 
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1)k (k + b-l)l In 
^ ' (k + b- P-n)\\k 



(-n- 
1 



(1 -x^-Pdx 



„b-l 



(-n- (3)1 Jo 



(l-x)-Pdx 



Putting ( |5.18[ ) into l |5.12[ l yields ( |5.13[ ) in this case too. 

Assume finally that n < A - 1 and /3 > b + 1. We write in ( |5. 1 2[ > that 

(fc + 6- 1)103 — ft + n-fc-l)! = C8 + n-l)! / x^^l - xf-^-^ dx . 



Then 

nA(iV-l) 



53 (fc + &-l)!(0-& + n-fc-l)!^J 



= (/3 + n- 1)! 



E 

, fc=0 



-6-1 



(l-x^+^dx 



= (0 + n - 1)1 J x'-Hl - xf-'- 1 dx = (6 - 1)109 - 6 - 1)! (/ ^! ^ ■ 

Putting ( f57T9] > into | (5TT2| | yields ( f57T3] >. 

By choosing /3 = 6 in Theorem 5.10 we derive that 

nA(JV-l) 



(5.19) 



□ 



We immediately obtain the following particular result which will be used in Theorem |6.1 
Corollary 5.11. The factorial pseudo-moments of (S£ — b) are given by 



4(#- & )J = 



\(-b) n ifO <n< A -1, 
I ifn > N. 



The above identity can be rewritten, if < n < A — 1, as 

2+ 



E 



■S+-b 



= (-iy 



n+b-l 
6-1 



(5.20) 



Moreover, since € {6,6+1, . . . ,6 + A-l}, it is clear that (S b + - 6) (S+-6-1) • • • (S,, 1 
6 — A + l) =0 which immediately entails that (S + — b) =0 for any n > A; th 



6 — A + l) =0 which immediately entails that (S^~ —b) = for any n > A; then 
E[(S^ - b) ] = for n > A as stated in Corollary 
By choosing /3 = in Theorem 



5.10 



5.11 



we plainly extract that E[(S£~) ] = for n € 
{1, . . . , A- 1}. Moreover, as previously, (S£) n = for any n > 6 + A; then E[(S^) J = 
for n > 6 + A. Actually, we can compute the factorial pseudo-moments of S + , E^S^") ] 



for n G {A, A + 1, . . . , 6 + A}. The formula of Theorem [5~To 
we provide another way for evaluating them. 



seems to be untractable, so 
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*[(#)»] 



Theorem 5.12. The factorial pseudo-moments of 5^ are given by 

(-if- 1 (Z'Ji) ( b + N - 1)« ifN<n<b + N- 1, 
ifl<n<N -lorn>b + N. 

Moreover, for n £ {1, . . . , N — 1}, the pseudo-moment of S£ of order n vanishes: 

E[(s+n=o 

and 



E 



TO 



N 



(-1) 



<-b) 



N ■ 



(6-1)! 

Proof. We focus on the case where N <n <b + N — 1. We have that 

JV-l 

E [(#) n ] - E p {^ + = 1 + b ^ e + b ^ 



=0 

b + N -1 
b 



N-l 



+ l)n-l 



iV-1 



The intermediate sum lying in the last displayed equality can be evaluated as follows: 
by observing that [l + b — l) n -i = {d 11 ^ 1 /dx n ~ 1 )(x l+b ~ 1 )\ x= i and appealing to Leibniz 
rule, we obtain that 

N-l 

£(-l)V + 6-l)n-l 



e=o 



N - 1 



dx 

n-l 



3i EH)' 



. 1=0 



N — 1 



€+6-1 



\JV-1 



dx r 
da;n-l-j 



x=l 



n-l 
- l 



= (-i^-^JV- 1) 
Consequently 

= (-1)"" 1 



(6-1)„_jv = (-1) 



AT-l 



(6-l)!(n-l)! 



6 + 7V- 1 
b 



(n- N)\(b + N -n-l)!' 
(6- l)!(n- 1)! 



(n-iV)!(6 + iV-n - 1)! 
(n -!)!(& + JV-1)! 



(n - N)\(N - 1)!(6 + N - n - 1)!' 

This is the result announced in Theorem |5 . 1 2 1 when N<n<b + N— 1. 

Next, concerning the pseudo-moments of 5^, we appeal to an elementary argument 
of linear algebra: the family (l,Xi,X 2 , ■ . ■ ,X n ) (recall that X k = X(X-1) ■ ■ ■ (X-k+1)) 
is a basis of the space of polynomials of degree not greater than n. So, X n can be 
written as a linear combination of Xi,X 2l X n . Then E[(S , h ^)™] can be written as a 
linear combination of the factorial pseudo-moments of of order between 1 and n. 
These latter vanish for n € {1, . . . , N - 1}. As a result, E[(S^) n ] = 0. 

The same argument ensures the equalities E[(5^ h ) JV ] = E 



(St 



N' 



, which is equal to 



(-l)"- 1 ^ + N- 1) N , and E[(S+ - b) N ] = E 
them yields the value of E (S£) N 



(St 



N 



(—6) which vanishes. Each of 



The proof of Theorem [5712] is completed. 



□ 
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6 Link with the high-order finite-difference operator 

Set A+f(i) = f(i + 1) - f(i) for any i e Z and (A+p = A+ o • • • o A+ for any j e 

j times 

N*. Set also (A + ) / = /. The quantity (A + p is the iterated forward finite-difference 
operator given by 



(a+)V60 = £(-i> 



fe=0 



Conversely f(i + k) can be expressed by means of (A+p f(i), < j < k, according as 

k 



3=0 



/(* + *) = E(JJ ( A+ ) J /«- 



(6.1) 



We have the following expression for any functional of the pseudo-random variable St . 
Theorem 6.1. We have, for any function f defined on {6, 6 + 1, . . . , b + N — 1}, that 



N-l 



e[/(#)] = £(-i) 



3=0 



j + b-1 

6-1 



(A+)V(6). 



(6.2) 



Proof By ( |6.1[ ), we see that 

'S+-6 



E[/(5+)]=E 



St - b 



N-l 



(A + y/(6) = E E 



S+-6 



which immediately yields ( |6.2| l thanks to ( |5.20[ l. 
Corollary 6.2. The generating function of is given by 



N-l 



E 



3=0 



3 + b-l 
6-1 



(i-C) j 



□ 



(6.3) 



Proof. Let us apply Theorem 6.1 to the function f(i) — ( l for which we plainly have 
(A + ) j f(i) = (-1) 3 C(1 ~ O 3 - This immediately yields joj. □ 

Remark 6.3. A direct computation with ( |5.10[ ) yields the alternative representation: 



E C b = b 



C° / x h - 1 {l-Cx) a - l dx 



N-l 



Of special interest is the case when the starting point of the pseudo-random walk is 
any point x € Z. By translating 6 into 6 — x and the function / into the shifted function 



/(. + x) in formula 1 6.2 1, we get that 



N-l 



E* [/(#)] = E[f(x + S+ J] = £(-1) 



j + 6 — x — 1 

i 



(A + )VW. 



Thus, we obtain the following result. 

Theorem 6.4. We have, for any function f defined on {6, 6 + 1, . . . , 6 + N — 1}, that 



N-l 



3=0 



(6.4) 
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ith P£ (x) = 1 and, for j e {1, . . . , N - 1}, 



1 



J ' k=0 

The < j < N — 1, are Newton interpolation polynomials. They are of degree not 
greater than (N — 1) and characterized, for any fee {0, . . . , N — 1}, by 

(A+)*P+(6) =*,•*. 

Proof. Coming back to the proof of Theorem |6.1| and appealing to Theorem |5.7| we 
write that 



P£,-(s) =E* 



S+ -b 



St„ + x - b 



N-l , v 
m=j ^ ' 



m\ {N —1 



(N- 



where, for any m G {0, . . . , N — 1}, 



K m (x) 



(6.5) 



K m (x) = ( Yl (b-x + k) \/{b-x + m) = Y[ (b-x + k). 

\k=0 ) o<fc<«-i 

The expression K m (x) defines a polynomial of the variable x of degree (N — 1), so PjJ" . 
is a polynomial of degree not greater that (N — 1). It is obvious that K m (b + 1) = for 
i£ {0,...,N -l}\{m}. On the other hand, K e (b + £) = (-l) e (N - ly./^J 1 ). By putting 
this into l|6.5[l, we get that 



KA b + £ ) = { j ) 1 {i<ty- 

Next, we obtain, for any k £ {0, . . . , N — 1}, that 

( A+ ) fcp M( & ) = E(-i) fc+£ (!) p ^( & + 4 = E(-!) fc+f A 

fc-3 



£=0 



fc - j 



The proof of Theorem|6.4|is finished. 



□ 



We complete this paragraph by stating a strong pseudo-Markov property related to 
time cr£ . 

Theorem 6.5. We have, for any function f defined on Z and any n e N, that 

N-l 

Mf&t+n)] = E p wW (A + y'E 6 [/(5 n )]. (6.6) 

i=0 



In 6.6), the operator (A+)- 7 acts on the variable b. 
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Proof. We denote by the pseudo-probability associated with the pseudo-expectation 
E x . Actually it represents the pseudo-probability related to the pseudo-random walk 
started at point x at time 0. We have, by independence of the Uj's, that 

Mf&t+n)] = E M%t=K,s+=e}f( S t + U k+ i + ■■■ + U k+n )] 

b<£<b + N -1 

= E ¥ ^ a i = k > s b = t}nf(t +u 1 + --- + u n )} 



fc,<€N: 
b<£<b + N-l 



b+N-1 



= E V X {S+ =£}E e [f(S n )}=E x [E s+ [f(S n )]]. 



i=b 



Hence, by setting g(x) = E x [f(S n )}, we have obtained that 

which proves l|6.6[l thanks to (|6.4[l. □ 



Example 6.6. Below, we display the form of ( 6.6 1 for the particular values 1, 2, 3 of N. 
• For N = 1, l |6.6| l reads 

which is of course trivial! This is the strong Markov property for the ordinary 
random walk. 



For N = 2, ( |6T6] > reads 

Mf(S at+n )] = E b [f(S n )} + (x-b) A+E b [f(S n )] 

= (b-x + l) E b [f(S n )} + (x-b) E b+1 {f(S n )}. 



• For N — 3, l |6.6| l reads 

Ex[/(^++n)] = E fc [/(5„)] + (x - b) A+E b [f(S n )} 

+ \(x-b)(x-b-l)(A+) 2 E b [f(S n )} 
= \ (x - b - l)(x - b - 2) E b [f(S n )} ~(x-b)(x-b-2) E b+1 [f(S n )} 
+ ^(x-b)(x-b-l)E b+2 [f(S n )\. 

7 Joint pseudo-distribution of (r 6 + ,X 6 + ) 

Below, we give an ad hoc definition for the convergence of a family of exit times. 

Definition 7.1. Let ((X^) t >o)s>o be a family of pseudo-processes which converges to- 
wards a pseudo-process (X t ) t > when e — > + in the sense of Definition 4.1 Let I be a 
subset ofR and set rf = inf{f > : Xf £ I}, Xf = X £ Tf and r z = inf{i > : X t £ I}, 
Xj = X Tj . 

We say that 



e— S-0+ 



if and only if 



VA > 0, V/i e E, E(c- A ^ +i ^l {T|<+oo} J — » E(e- A ^ +i ^l {Tj<+oo} ^ 
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We say that 
if and only if 



xi —> x, 



V/i e R, E(e i ^l {r|<+oo} ) — » E(e f ^l 



1 {t j <+do} 



As in Section[4] we choose for the family ((X t £ ) t > ) £ >o the pseudo-processes defined, 
for any e > 0, by 

and for the pseudo-process (X t )t>o the pseudo-Brownian motion. For I, we choose the 
interval (-00,6) so that rf = r b £+ , Xf = X e b + and Tj = r+, X 7 = X+. Set & £ = \b/e\ 

where [ . ] is the usual ceiling function. We have = e 2N ' and X^ + = eS^ . Recall 

• 2 j — 1 

the setting tpj = — ic 1 TJ^ 7r / 1 < j < N. 

Theorem 7.2. Assume that c < 1/2 2N ~ 1 . The following convergence holds: 



e— >0+ 



where, for any A > and any /x e 



{t(,<+oo} 



= C 



fc = l 1<3<JV ft ft l<j<JV V V 



Proof. We alrea dy pointed out that the assumption c < 1/2 27V 1 entails that = 1. 
Therefore, 1 5.7 1 holds for z = c~ Ae " < l/M^ = 1, i.e., for A > 0. So, by (5.7 1, we have, 
for A > 0, that 



K+<+oo} 



-A £ 21v o-+ +i/ieS+ 

2JV 



{<T+ < + Oo} 



£L fc (e-^,e-) Ufc (e-^) ; 



fe=i 



Recall that we previously set Uj (A, e) = Uj (e Ae " ) . Thanks to asymptotics 1 4.5 1, we get 
that 

u k (X,e)-u 3 (X,e) ~ - cp k y<fx/c'e, 1 — ttj(A, e)e l ^ e ~ (^/X^ - i/i) £ . 

e— >0+ e— >0+ 



Thus, 



►0+ 



lim L fe (e- AE2 ",e^)= JJ 



TT ft TT h_i^j£ 

^4, (ft ~ <^)Wc w - ft M» v 2 v^A 



Finally, we can easily conclude with the help of the elementary limits 

lim e^ £b ' =e^ b , lim u k (\,e) b ? = e^'v^ 6 . 

£-►0+ e->-0+ 



□ 



Theorem 7.3. The following convergence holds: 
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where, for any ^iSl, 



j=o J 

This is the Fourier transform of the pseudo-random variable X^. Moreover, 

P{r+ < +00} = 1. 

Proof. By ( |6.2[ ), we have that 

E ^ + l K+<+ ^ } )=E(e^<l { , <+oo} 



We can easily conclude by using the elementary asymptotics 

lime^^e'" 4 , fj+b £ -l\ _ V ^_ ifl ey „ (_ i//e )J. 

□ 

Corollary 7.4. The pseudo-distribution of X^ is given by 

This formula should be understood as follows: for any (N — l)-times differentiable 
function /, by omitting the condition < +00, 



E[/W)] = E(- 1 ) i ir /W(6) - 

3=0 ■'' 

We retrieve a result of 1 1 ] and, in the case N = 2, a pioneering result of 111 711 
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Part III — First exit time from a bounded interval 



8 On the pseudo-distribution of (a a b, S a b) 

Let a, b be two integers such that a < < b and let £ = {a, a — 1, . . . , a — N + 1} U 
{b, b + 1, . . . , b + N — 1}. In this section, we explicitly compute the generating function 
of (a ab , S a b). Set, for I e £, 

H ab/ (z) = e(^i {So6=Acto6<+(X)} ) = Yl p K*< = k > s -b = e}* k - 

fc6N 



We are able to provide an explicit expression of H ab ^(z). As in Section[5] due to 12.14 1, 
we have the following a priori estimate: |P{c a & = k, S ab — £}\ < |P{<Si € (a, b), . . . , Sk-i € 
(a,b),Sk — £}\ < M*. As a byproduct, the power series defining H a b/(z) absolutely 
converges for \z\ < 1/M\. 



8.1 Joint pseudo-distribution of (a ab , S ab ) 

Theorem 8.1. The pseudo-distribution of (<T ab ,S ab ) is characterized by the identity, 
valid for any z e (0, 1), 



E(z' Tab l 



where 



1 Ul 



Wi{ui(z), . . - ,u 2 n(z)) 
W(ui(z), . . .,u 2 n{z)) 



(8.1) 



W( Ul 



,U 2 N 



N-1 b-a+N-1 b-a+2N-2 
1 Ui ■ ■ ■ i*i 



b-a+N-1 

Urn AT 1*0 



1 U 2 N ■ ■ ■ U 2N U 2N ■ ■ ■ U 2N 

and, if a — N + 1 < £ < a, Wg(ui, . . . , U2n) is the determinant 

1 Ul 



b-a+2N-2 



e+N-a-2 „N-a-l 



II, 



l+N-a 



N-1 b-a+A»-l 



f + jV- a _2 N-a-1 l+N-a 



1 U2N ■ ■ ■ <^2N U 2N "2JV • • • "27V 

if b < I < b + N — 1, We(ui, . . . , u 2 n) is the determinant 

1 Ml 



U 



„N-l b-a+N-1 

U"3 AT t*Q 



*2N 



N-1 ^b-a+N-1 



1 "27V 



u 



b-a+N-1 
2N ° 



E+N-a-2 „,N-a-l „,£+N-a 



U 



l 2N 



P+N-a-2 N-a-1 l+N-a 
2N a 2N a 2N 



b-a+2N-2 



6-a+2JV-2 
l 2N 



,b-a+2N-2 



b-a+2N-2 
l 2N 



Proof. Pick an integer k such that k < a or k > b. If S n = k, then an exit of the 
interval (a, b) occurs before time n: a a b < n. This remark and the independence of the 
increments of the pseudo-random walk entail that 



¥{S n = k}= ¥{S n = k,a ab <n} = J2J2 = fc ' a ">> = S - b = ^ 

j=a ies 

n 

0=0 ies 



(8.2) 
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Thanks to the absolute convergence of the series defining Gk(z) and H ab ^(z) for z <= 
(0, 1) and \z\< 1/K'h respectively, we can apply the generating function to eguality ( 8.2 i. 
We get for z € (0, 1/Mi), that 

G k (z) = J2 G k-e(z)H ab j(z). 



Using expression (3.5 i of Gk, namely Gk(z) — YljLi a j( z ) we 9 e ^ f° r k > b 

N - 1 (recall that vj{z) = l/uj(z)), that 



J=l \t££ I 



and, for k < a - N + 1, that 



(8.3) 



(8.4) 



When limiting the range of k to the set {b+N, 6+JV+l, . • • , 6+2iV-l} in ( (Oj l and to the 
set {a-2iV+l, a-27V+2, . . . , a~iV} in ( |8.4| l, we see that \8.3\ and l |8.4[ l are homogeneous 
Vandermonde systems whose solutions are trivial, that is, the terms within parentheses 
in (|8.3p and l|8.4[l vanish. Thus, we get the two systems below: 



^H aM (z)u,(z) £ = l, l<j<N. 



£<££ 



and 



ee£ 



It will be convenient to relabel the Uj(z)'s and vj(z)'s, 1 < j < N, as «j(.z) = U/+jv(z) and 
= itj+jv(z); note that Vj(z) = l/uj(z) for any j e {1, . . . , 2^} and {wi(z), . . . , u 2 n(z)} 
= {vi(z), . . . ,v 2 n(z)}. By using the relabeling Uj(z),Vj(z), 1 < j < 2N, we obtain the 
two eguivalent following systems of 2N equations and 2N unknowns, U\(z), . .. ,u 2 n(z) 
for the first one, ui(z), . . . ,V2jv(-?) for the second one: 



2fr oW («)« i (2) / = i, i<j<2iv, 



£<££ 



and 



^ Hab^Vjizf = 1, 1 < i < 2JV. 



(8.5) 



(8.6) 



Systems (8.5 i and \Q.Q\ are "lacunary" Vandermonde systems (some powers of Uj(z) are 

1 < j < 2N. 



missing). For instance, let us rewrite system (|8.5p as 



(8.7) 



l££ 



Cramer's formulae immediately yield (8.1 1 at least for z e (0, 1/Mi). By analyticity of 
the Uj's on (0, 1), it is easily seen that (8.1 1 holds true for z e (0, 1). Systems (8.5l and 
( |8.6| l will be used in Lemma [83] □ 

A method for computing the determinants exhibited in Theorem |8.1| and solving 
system i 8.7[ l is proposed in Appendix |A.1| In particular, we can deduce from Proposi- 
tion A.3 an alternative representation of E(z <Tab lis a b=i,^ab<+oo}) which can be seen as 
the analogous of ( |5.1[ ). Set s (z) = 1 and, for k,£ G {1, . . ., 2N}, 

si(z) = ^ Ui^z) ■ ■ -u H (z), p k (z) = Yl i u k(z) - Ui(z)], 

Ki 1 <-<i t <2N l<i<2JV 
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s k,o{ z ) = h s k,m(z) — for any integer m such that m < — 1 or m > 2N and, for 
m € {!,..., 27V- 1}, 



l<i 1 < <i m <2N 
i 1 ,...,i ml tk 



Set also 



W(z) 



W ke (z) = 



s N (z) s N -i(z) 
s N+1 (z) s N (z) 

SN+b-a-2(z) SN + b- a -3(z) 

Sk,N{z) S k ,N-l(z) 
Sk,N+l(z) Sk,N{z) 

Sk,N+b-a-2{z) Sk.N+b-a-3{ z ) 

Sk,N+b-e-i(z) su.N+b-e.-2{z) 



SN-b+a+2(z) 
SN-b+a+l(z) 

s N (z) 

Sk,N-b+a+l(z) 
Sk,N-b+a{z) 

Sk,N-l(z) 
Sk,N+a-e(z) 



Then, applying Proposition A. 3 with the choices p = r = N and q = b — a — 1 leads, for 
any I G £ , to 



E(z^l {Sab= ^ ab<+oo} ) 



(-1) 



l+N-a-l 2N u 



\- W M {z) 



N-a-1 



(8.8) 



w (*) k~i Pk ^ 

The double generating function defined by 

E(z^( s ^l {aab<+oo} ) = ^E(^l {Sa6=Va6<+oo} ) C l 

admits an interesting representation by means of interpolation polynomials that we 
display in the following theorem. 

Theorem 8.2. The double generating function of (a abl S ab ) is given, for z e (0, 1/iWi) 
and C e C, by 



2N 



E(z^( s «n {aab<+oo} ) = £*(*, (v k (z)0 

fc=l 



a-jV+1 



(8.9) 



where 



l k (z,C) = Pk(z,() II 7 \ ' ke{l,...,2N} 

are interpolation polynomials with respect to the variable C satisfying L k (z, uj(z)) = Sjk 
and Pk{z, C) are some polynomials with respect to the variable C of degree (b — a — 1). 



Proof. By (|8.1[), we have that 



E(z""><; s «n {Sab < atlTab<+oo} )= H ab,i{^= E 



We(u!(z), . . .,u 2 n(z)) h 



e =a -N+l 



=a^N +1 W(u 1 (z),...,U 2N (z)) 



In order to simplify the text, we omit the variable z. We expand the determinant 

We(ui, . . . , U2n), a — N + 1 < £ < a with respect to its {I + N — a)th column: 



2N 



W e (ui, . . . ,u 2N ) = X] M f ° 1 W ke (ui,...,U2 N ) 



k=i 
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where W u {u u 

1 Ul 



, u 2n). 1 < k < 2N, is the determinant 



l+N-a-2 



u\ 



l+N-a 



1 


Uk-1 ■ 


e+N-a-2 
"fe-1 





it 


1 


u k 


l+N-a-2 
u k 


1 


u 


1 


Uk+1 ■ 


l+N-a-2 
"fc+1 





u 



e+N-a 
k-l 



i+N-a 
fc+1 



1 u 2N 



e+N-a-2 
l 2N 



u. 



2N 



a 




6-o+JV 
"l 


a 




"fe-1 


a 




U fc 


a 




"fc+1 


a 


U 2JV 


,.b-a+/V 



which plainly coincides with 



1 Ul 

1 Ufe_l .. 

.. 

1 u k+ i . . 

1 U 2 N ■ ■ 



U 



e+N-a-2 o e+N-a-l 0l i+N-a 



U 



u 



e+N-a-2 ^e+N-a-1 n ,i+N-a 



k-l 



u 



k-l 







1 



fe-i 




u 



e+N-a-2 o e+N-a-l 0l i+N-a 



fc+1 



U 



fc+1 



fc+1 



u 



e+N-a-2 ^e+N-a-1 n ,i+N-a 



2N 



u. 



2N 



Therefore, we obtain that 



27V 



' 2N 



N-l n b-a+N-l 
1 Ui 



,JV-1 h-a+N-1 
L k-1 u k-l 



N-l h-a+N-1 
l k+l U k+1 



N-l h-a+N-1 
HN u 2N 



b-a+2N-2 



b-a+2N-2 
l k-l 

b-a+2N-2 
l k 

b-a+2N-2 
l k+l 



b-a+2N-2 
l 2N 



,b-a+2N-2 



,6-o+2/V-2 
*fc-l 



b-a+2N-2 
l k+l 



b-a+2N-2 
l 2N 



Y, W e ( Ul ,...,u 2N )C e = J2 [Y u k~ a ~ lw ^i,...,U2N))C 

i=a-N+l e=a-N+l \fc=l / 

2N / a 

= c -n+ij2( Y w kl (ui,..., U2N )c 

fc=l \e=a-N+l 



i+N-a-1 \ y a-N+l 



Next, we can see that the foregoing sum within parentheses is the expansion of the 
following determinant with respect to its fcth row (by putting back the variable z): 



1 ui(z) ... u\ 



b-a+N-l 



(z) 



b-a+2N-2 



(z) 



Uk-l(z) 

c 

u k +i(z) 



■ <-i'{z) <__\ +N -\z 



u 



N-l 



c 

N-l 
l k+l 



(z) U 





b-o+JV-l 
fc+1 



b-a+2N-2( \ 
k-l \ Z > 





(z) 



u 



b-a+2N-2 
fc+1 



1 u 2N (z) ... u^{z) ut N a+N -\z) 
As a result, by setting 

D(z) = W{ui{z), . . . ,u 2N {z)) 

1 ui{z) ... u?-\z) u\- a+N - 1 {z) 

1 u 2N (z) ... u^\z) u h 2N a+N - l (z) 



b-a+2N-2 
l 2N 



(z) 



(z) 



b-a+2N-2 



,b-a+2N-2 



J '2N 



(z) 



(z) 
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we obtain that 



E(^C^^<«,<^<+~}) = J2 D k^0(vk(z)Q a - N+1 - (8-10) 



Similarly, we could check that 

E(z^C Sab Ms ab >b,a ab <+c 

where D k (z, C) is the determinant 



>) = J2 D i( z >0Mz)cr- N+1 (8.H) 



fc=i 



1 ui(z) ... uf _1 (z) u b 1 - a+N - 1 {z 



... ^b-a+W-l 

1 u k+l {z) ... u™-t{z) ul-^iz 



1 u 2N (z) ... u^\z) u h 2N a+N - l (z 



By adding < |8.10[ l and l |8.11[ l and setting 

D k (z,C)=D+(z,0 + D^(z,C) = W(u 1 (z),...,u k . 1 (z 
we obtain that 

2N 



b-a+2N-2 



b-a+2N-2 
1 k-l 



(z) 



"b-a+2N-2 



b-a+2N-2 
1 k+l 



(z) 



4 N a+2N - 2 (z) 



X,u k+1 {z), . . .,u 2N (z)), 



E(z^{ s °n {aab<+oo} ) = J—^D k (z,Q (v k (z)C) 

^ ' k=l 



a-N+1 



(8.12) 



We observe that the polynomials L k (z, () = D k {z, ()/D(z) with respect to the variable ( 
are of degree b — a + 2N — 2 and satisfy the egualities L k (z, Uj (z)) — 5j k for all j £ £ and 
k e {1, . . . , 2iV}. Hence they can be expressed by means of the elementary Lagrange 
polynomials as displayed in Theorem|8.2| □ 



Example 8.3. For N = 2, ( [8T9] > reads 

E(z^( s °n {aab<+oo} ) =U{z,0 Mz)C) a ~ J + U(z,0 (V2(z)0 



a—1 I t" / >\ / / \ >\a— 1 



+ L 3 (z, C) (v 3 (z)C) a - L + L 4 (z, C) (v 4 (z)CT 



where 



- . = W(C,u 2 (z) 1 u 3 (z),u 4 {z)) ~ _ W(u 1 (z)X,u 3 (z),u 4 (z)) 

W{ui{z),u 2 {z),u 3 {z),UA c (z)y 

r i r \ - W ( u l( z )> u 2(z),(,M z )) f / ,s 

^{Z, (J — —r — — — 7-rr, L, 4 (Z,(,) 

W[ui{z),u 2 {z),u 3 {z),u A (z)) 



W(ui(z), u 2 (z), u 3 (z), Ui(z)) ' 

W(u 1 (z),u 2 (z),u 3 (z)X) 
W(ui(z), u 2 (z), u 3 (z),u A {z)) ' 



All the polynomials L k (z, (), 1 < k < 4, have the form A kta ^ 1 (z)+A k , a {z)C,+A k ^ b (z)(; b - a+1 
+A k , b+1 (z)( b - a +\ 
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Remark 8.4. By expanding the determinant D k (z,() with respect to its fcth raw, we 
obtain an expansion for the polynomial L(z, () as a linear combination of 1, (, ... , C^ -1 , 

£b-a+N-l^b-a+N^ ^ ( *b-a+2N-2_ that ^ an expansiori f t he form 

AT-1 b-a+2N-2 

1=0 e=b~a+N-l le£ 

Hence, 

2N ( \ (IN \ 

E^-c^-w+co}) = E E^m*)^ Mz) a ~ N+1 = E E^hw^ 1 c £ 

k=i \e<££ ) ie£ \k=i ) 

from which we extract, for any I e E, that 

2JV 

E(z^t {Sab ^ aab<+oc} ) =Y,Au{z)u h {z) N - a - 1 . 

k=l 

Actually the foregoing sum comes from the guotient We(ui(z), . . . , U2n(z))/W(ui(z), . . . , 



u 2 n(z)) given by 1 8.1 1 by expanding the determinant Wi(ui(z), . . . , u 2 n(z)) with respect 
to the [l + N — a)th column or (£ + N — b + l)th column according as a — TV + 1 < I < a 
orb<£<b + N-l. 

8.2 Pseudo-distribution of 5 a6 

In order to derive the pseudo-distribution of 5 ab which is characterized by the num- 



bers H a bj{l), I e £, we solve the systems obtained by taking the limit in 1 8.12 1 as 

z -> 1". 

Lemma 8.5. The following identities hold: for N < k < 2N — 1, 



t=b 



^ fb + N -i-e\ TT _ (b + N-i\ tOAAS 

E h Rm(i)= fc ■ (8 - 14) 



£=n-JV+l 



Proof. By 1 3.3 1, we have the expansion = 1 + £j{z) where £j{z) = 0( 2 \/l — z) 

for any j e {1, . . . , N}. Actually such asymptotics holds true for any j e {1, . . . , 27V} 
because of the eguality uj = l/uj-N for j € {N + 1, . . . , 2iV}. We put this into systems 
(|8.5p and l|8.6[l. For doing this, it is convenient to rewrite these latter as 



YHabAz)u J (z) e+N - a ~ 1 = Uj (z) N - a -\ l<j< 2N, 
ee£ 

Y H ab,l{ z ) u i{ z ) h+N ~ 1 ~ l = ^{z) b+N -\ 1 < j < 2N. 



ie£ 

We obtain that 



J2^+^(z)Y +N - a - 1 H ab Az) = (l+s ] (z)) N - a -\ l<j<2N, (8.15) 



l<££ 



J2(l + £j(z)) b+N - 1 - e H abA (z) = (l + e j (z)) b+N -\ l<j<2N. (8.16) 
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System ( |8.15p writes 

E 



b-a+2N-2 

E 

fc=0 



t>k + a-N + l 



e + N-a-1 
k 



H a bA Z ) £ A Z ) h 



N-a-1 

E 

k=0 



N-a-1 
k 



ej{z) k - 



(8.17) 



Set 



M k (z) 



e<=£: 
£>k + a-N + l 



b-a+2N-2 



N-a-1 

k 



E 

+ 

J2 M k (z)ej(z) 



H a b,e( z ) - 



N-a-1 
k 



k=2N 



Then, equality l |8.17| l reads 

2N-1 

M k (z) £j (z) k = Rj(z), 1 < j < 2N. 



k=0 



This is a Vandermonde system which can be solved as in the proof of Lemma [5^6] upon 
changing N into 2N. We can check that lim 2 _ > . 1 - M k {z) = 0, which entails that 



E 



£>k + a-N + l 



£ + N -a- 1 
k 



H a b : e(l) 



Similarly using ( |8. 1 6p , we can prove that 



E 



£<b + N -1-k 



b + N -£-1 
k 



H a b,e(l) — 



N-a-1 
k 



N + b- 1 
k 



, 0<k<2N-l. 



, < j < 2N - 1. 



(8.18) 



(8.19) 



Actually we shall only use \Q.1Q\ and l |8.19| l restricted to j € {N, . . . , 2N - 1} which 
immediately yields system ( |8.13 1-( |8.14[ ). □ 

Now, we state one of the most important result of this work. We solve the famous 
problem of the "gambler's ruin" in the context of the pseudo-random walk. 

Theorem 8.6. The pseudo-distribution of S a b is given, for I G {0, 1, . . . , N — 1}, by 



¥{S ab = a-£,a ab < +c»} = (-1)' 



KN fN - 1\ ,/l + b-a + N-l 



I- a 



N 



V{Sab = b + £, a ab < +00} = (-1)' 



KN / N — 1\ Jl+b-a + N-1 



l + b 



N 



where 



K 



N-a- l\(N + b- 1\ (-1) 



N 



N 



AH 



• a(a - 1) . . . (a - N + 1) 6(6 + 1) . . . (6 + N - 1). 



Moreover, ¥{a ab < +00} = 1 and 



¥{<t+ < a'} = KN 2 / / u- a -\l - u) N - 1 v»- 1 (l - v) N ^ dudv 



v- 



Oa <<y b } = KN 2 // u- a ^ i (l-u) N - i v b - i (l-vy^- i dudv 1 



v- 



From pseudo-random walk to pseudo-Brownian motion 



41 



where 



V + = {{u,v) eR 2 :0<»<(i<l}, V~ = {(u,v) e R 2 : < u < v < 1}. 

Proof. We have to solve system l |8.13[ l-l [8.14[ l. For ( |8.13[ l for instance, the principal 
matrix and the right-hand side matrix are 



+ b-a + N - 1 
k + N 



and 



J o<k,e<N-i 

and the matrix form of the solution is given by 

i + b-a + N - 1\ ] 1 

k + N J \o<k.i<N-i 



N-a-1 
k + N 



0<k<N-l 



N-a-1 
k + N 



0<k<N~l 



The computation of this product being quite fastidious, we postponed it to Appendix |A.3| 
The result is given by Theorem |A.8| 



(-1)' 



KN (N - 1\ ,(l + b - a + N - 1 



e + b\ i 



N 



0<1<N-1 

The entries of this matrix provide the pseudo-probabilities P{S a b = b + l,a a b < +oo}, 
< i < N — 1, which are exhibited in Theorem 8.6 The analogous formula for P{S a b — 
a — I, a a b < +00} holds true in the same way. 

Next, by observing that a ab = rnm(a~,a£) and that {a ab < +00} = {a~ < +00} U 
{cr^ < +00}, we have that 

P K + < °a } = F K + < °a'°ab < +«^} = V{S ab > b,O ab < +CX)} (8.20) 
N-l 



W{S ab = b + £,a ab < +00} 



N-l 



N {N-l\,{i + b-a + N-l 



N 



(8.21) 



Noticing that l/(i + b) = / X y 1 ^ 1 dy and l/(' +b " a w +JV " 1 ) = Nj^ - x )N-i dXj 

we get that 



N-l 



EC- 1 ) 



N fN - 1\ , (i + b - a + N - 1 



e=o 



£ + b\ i 



N 



1 ,1 /N-l 



N 1 



N 1 



0J0 
1 i-i 



EC- 1 ) 



x b - a - 1 {\ - x) N - L y b - L {\ - xy) N - L dxdy 



N-l 



(xyf 1 x b - a - 1 (l-x) N - 1 y b - 1 dxdy 



\N-l„.b-l 



N-l 



JO 



The computations can be pursued by performing the change of variables (u, v) = (x, xy) 
in the above integral: 



N-l 



£(-*) 



N fN - 1\ , (I + b - a + N - 1 



e=o 



£ + b\ i 



N v 



N 



u- a - l {l - u) N - l v b ~\l - v) N ~ l dudv 



(8.22) 



v+ 
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Putting ( 8 . 2 2 i into (8.211 yields the expression of P{ a£ < a a } displayed in Theorem 
The similar expression for P{a~ < a b } holds true. Finally, 

P{cr Q b < +00} = P{cr- < (T+,<7 ab < +00} + P{cr+ < a~, a ab < +00} 



8.6 



H< <a+} + F{<j+ <a-} 



KN 



u 



^l-ttr-v-^i-w) 



N-l 



du dv. 



0J0 



The foregoing integral is quite elementary 

-1 r i 

u 



-a-l 



(l-«) 



7V-1 



du 



B(N,-a)B(N, b) 



v^il-vf-'dv 

1 1 



which entails that P{(7 ab < +00} = 1. □ 

In the sequel, when considering S ab , we shall omit the condition a ab < +00. 
Example 8.7. Let us have a look on the particular values 1, 2, 3 of N. 

• Case N = 1. In this case S ab £ {a, b} and 

¥{S ab = a} = P{a~ < a+} = ¥{S ab = b} = ¥{a+ < a~} = 

We retrieve one of the most well-know and important result for the ordinary ran- 
dom walk: this is the famous problem of the gambler's ruin! 

• Case N = 2. In this case S ab e {a — 1, a, b, b + 1} and 

ab{b+l) 



¥{S ab = a - 1} = 



(6- a + !)(&- a + 2)' 



nSab = a} = - 



(fl- 1)6(6+1) 
(6-a)(6-a + l)' 

a(a — 1)6 

(6-a)(6-a + l)' " L "" u ~ ' ~> (6-0 + l)(6-a + 2)' 
6(6+ 1)(6- 3a + 2) 



p{^=6}= " (a -y b+i L p{^=6+i } 



PR" < *6 } = 
PR < <Ta } = 



'' 1 ~ (6-a)(6-a + l)(6-a + 2)' 
a(a- l)(36-o + 2) 
(b-a)(b-a + l)(b-a + 2)' 

Case N = 3. In this case S ab e {a — 2, a — 1, a, 6, 6 + 1, 6 + 2} and 

a(a- 1)6(6 + !)(& + 2) 



P{S ab = a - 2} 
P{S ab = a - 1} 
P{5 a6 = a} 

P{S a6 = &} 
P{S ab = 6+l} 
¥{S ab = 6 + 2} 



(6 - a + 2) (6 - a + 3) (6 - a + 4) ' 
o(a - 2)6(6 +1)(6 + 2) 

_ (6 - a + 1)(6 - a + 2) (6 - a + 3) ' 

(a - l)(a-2)6(6+l)(6 + 2) 

(6-a)(6-a+l)(6-a + 2)' 
0(0- l)(a-2)(6 + 1)(6 + 2) 

~(6-a)(6-a + l)(6-a + 2)' 
a(a- l)(a-2)6(6 + 2) 

(6-a + l)(6-a + 2)(6-a + 3)' 
a(a- l)(o - 2)6(6 + 1) 

~ (6 - a + 2)(6 - a + 3)(6 - a + 4) ' 
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+ 6(6+ l)0 + 2)(10a 2 -bab + b 2 -25a + 7b + 12) 
Wa <° b i- (fe _ a )(6 _ o + 1)(6 - o + 2)(6 - a + 3)(6 - a + 4) ' 
a(a - l)(a - 2)(a 2 - 5a& + 10& 2 - 7a + 25k + 12) 
' (b - a)(b - a + 1)(6 - a + 2) (6 - a + 3)(6 - a + 4) ' 



< °a } 



8.3 Pseudo-moments of S„ 



Let us recall the notation we previously introduced in Section 5.3 (i) n = j(i — l)(i — 
2) • • • (i — n + 1) for any i e Z and any n e N* and (i)o = 1, as well as the conventions 
1/i! = for any negative integer i and J2i=i = if i > j. 

In this section, we compute several functionals related to the pseudo-moments of 
Sab- Namely, we provide formulae for E[S a b(S a b ~ &)n-i] (Theorem 8.9 1, E[(S a (,)™] (Corol- 
lary 8.10i and E[(5 a& — b) n ] (Theorem |8. 11 1. This schedule may seem surprising; actu- 
ally, we have been able to carry out the calculations by following this chronology. 

Putting the identities l/( e+b ' a ^ N ~ 1 ) = N /„* ^-"^(l - x)*- 1 dx and l/(£ + b) = 
Jo y 1 ^ -1 dy into the eguality 



E 



e=o 



i + b V t 



£ + b-a + N - 1 
TV 



we immediately get the following integral representations for E[f(S ab )l{ Sab > b }], and 
the analogous ones hold true for E[f(S a b)lfs ab < a }] ■ 
Theorem 8.8. For any function f defined on £, 



E [/(£,&)! 



{S ab >b}\ 



= KN 



N-l 



KN 2 



fl 

10 JO 



N-l\ f(£ + b) t 
£ J £+b x 



f{£ + b){xyf 



EC- 1 

'N-l 

EM 

E[f(S ab )l {Sab <a}] 
JtJ l£=0 



„b-a-l 



\N-1 



dx 



(8.23) 



N-l 

£ 



(1-xY 

b - a - 1 (l-x) N - 1 y b - 1 dxdy, (8.24) 



„b-a-l 



(1-X) 



N-l 



dx 



= KN 1 



fl 

10 JO 



N-l 



E(-^ 



1=0 



N-l 



f(a-£)(x V y 



.b—a-l 



(1-X 



N-l„.-a-l 



(8.25) 



dxdy. (8.26) 



In view of ( |8.23| l and ( |8.25[ l and in order to compute the pseudo-moments of S a b, it 
is convenient to introduce the function /„ defined by f„ (i) = i(i — &) ra -i for any integers 
i and n such that n > 1. In particular, fx(i) = i. We immediately see that, by choosing 
/ = fi in Theorem |8.8| quantities ( |8.23[ l and < |8. 2 5| ) are opposite. As a by product, 
E[5 a b] = 0. More generally, we have the results below. 

Theorem 8.9. For any positive integer n, 



E[S a b(S a b - b) n -il{s ab >b}] 

\-ir- KN 



N\ 



2N -n (N-n)\ 



2N + b - a - 2 
27V -n 







if 1 <n<N, 
ifn >N + 1, 



(8.27) 



^•[S a b(S a b — b) n -l^{S ab <a}] 
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. KN N\ j2N + b-a-2 

(~ ) OAT ~ I AT Zm / 



/ n + 6 — a — 2 
n-2iV 



27V -n (iV-n)! ' \ 2iV - n 

(— l)- w+ " -x jRTJVJV! (n - N - 1)1 ( 

In particular, for any tie {1, . . . , 2N — 1}, 

E[S ab (S ab - &)„_!] = 0. 

Proof. By ( [8T23] >, we get that 



if 1 < n < N, 

ifN+l<n<2N-l, (8.28) 
if n > 2N. 



E[f n (S ab )l {Sab > b} ] — KN 2 



Mn-i a; 



t+b-a-l 



By noticing that (£)„_! = for £ € {0, . . . , n - 2} and ("j 1 ) (£) n _ a - (AT - l)„_ x 

for £ > n — 1, we obtain that 



N-l 



1=0 



N-l 



N-l 



7V-n 
i-n + 1 



N-r 



(N - l)^!^^ a;"" 1 £ (-1)^-1 

-1) 



1=0 



N -n 
I 



n - l( - N ^ x n - 1 (l-x) N - n if I <n<N, 



(N - n)\ 



if n > N+ 1. 



Hence, if 1< n < JV, 



= j KNNl x (n + 6-a-2)!(2JV-n-l)l 



(N — n)\ 



(2N + b-a-2)\ 



and we arrive at (8.27 1. Moreover, if n > N + 1 and 5 a b > b, we have S a 6 € {6, & + 1, 
. . . , 6 + N — 1}. Then, it is clear that f n (S ab ) = and 1 8.27 1 still holds in this case. 
On the other hand, by Theorem |8.6| we get that 

N-l 



E[fn(Sab)l {Sab <a}] = £ P{^a 6 = - £}/«(« - *) 



JV-1 



€=0 



iV 



iV-l 



= (~i) n KNNi j^^y 



1=0 



N-l\ {£ + b-a + n-2)\ 



£ J (e + b-a + N- 1)!' 



For n < iV, we can write that 

{i + b- a + n - 2)! 



+ 6 - a + N - 1)! (N - n)\ (£ + b - a + N - 1)( 



1 



(N-n)lJ 



A+n+b-a-2 



l+b-a+N-2\ 
N-n ) 



(l-x) N ' n dx 
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Then, 

= ^ x"+"- Q - 2 (l - a;)^-"- 1 dx 

KNN) (n + b-a-2)[(2N -n-l)\ 
~ ( ~ ' (N - n)\ X (2JV + &-a-2)! 

which proves l |8.28[ l. For n > N + 1, we write instead that 

(£ + b-a + n-2)\ _ d™-^ 1 
(£ + b-a + N-l)\ ~ da;"-^- 1 

Therefore, 

i-N-l 



x e+n+b-a-2 (1 _ x - ) N-n dx 



x l+n+b-a-2\ 



x=l 



in-N-1 "-^ /AT_1\ 

E[US ab )l {Sa ^ a} ]=(-irKNNl^^i^' i y{ £ ) 



JV-l 

f+n+b-a-2 



.T 



=0 

dx r ' 



(—l) n KNN\ (x«+ b — 2 (l x)^) 



x=l 



\i N + 1 < n < 2N — 1, the above derivative vanishes since 1 is a root of multiplicity 
N — 1 of the polynomial x e+n+b ~ a ~ 2 (l — x) 1 * -1 . Finally for n > 2N, we appeal to Leibniz 
rule for evaluating the derivative of interest: 



" X 1 k fn-N-l\ (n + b-a-2y 



x = l 



k=0 



(_l)"-i(JV_l)! (n + 6-a-2)! fn-AT-1 



(-l^-^n-iV- 1)1 



(2iV + 6-a-2)! V N-l 
n + 6 — a — 2^ 
n- 27V 



This proves < |8. 28p in this case. □ 
Corollary 8.10. Forn e {1, . . . ,2N— 1}, the pseudo-moment of S a b ofordern vanishes: 

n(s ab ) n ] = o. 

Moreover, 

E[(S ab ) 2N ] = -a(a — 1) • • • (a — -/V + 1)6(6 + 1) • • • (b + iV - 1). (8.29) 

Proof. As in the proof of Theorem [5712] we appeal to the following argument: the poly- 
nomial X n is a linear combination of fi(X), . . . , f n (X). Then E^S^;,)"] can be written 
as a linear combination of E[fi(S a b)}, ■ ■ ■ , E[f n (S ab )] which vanish when 1 < n < 2N — 1. 
Thus, E[(S' a& ) n ] = 0. The same argument entails that 

E[(S ab ) 2N ] =E[f 2N (S ab )] = (-If-'KM 2 

which proves l |8.29[ l. □ 
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In the theorem below, we provide an integral representation for certain factorial 
pseudo-moments of (S ab — a) and (S ab — b) which will be used in next section. 
Theorem 8.11. For any integer n e {0, . . . , N — 1}, 

E[(S ab - b) n l {Sab > b} ] = (-l) n KN 2 (N-l) n 

- a -\l - V+^l - v) N - n ~ x dudv, (8.30) 



E[(a - S ab ) n l {Sab < a }] = (~l) n KN 2 (N - 1)„ 

u n - a -\l - u)*-"- V-^l - vf- 1 dudv. 

v- 

The above identities can be rewritten as 



(8.31) 



E 



E 



Sab - b 

n 



a — S ao 
n 



{S ab >b} 



= {-l) n KN 



v+ 



{S ab <a} 



{-l) n KN 



N-r 

n 

- a -\l - u) N - 1 v n+b - l (l - v) N ~ n ^ dudv, (8.32) 

n- r 

n 

n ~ a -\l - uf-^v^il - vf- 1 dudv. (8.33) 



v- 



Proof. By ( |8.24p , we have that 

E[(S ab -b) n l {Sab > b} ] 



KN 



N-l 



EM) 



1=0 



N - 1 



(iUxyY 



„6-a-l. 



l-x) N - l y b - l dxdy. 



The sum lying in the above integral can be easily calculated: 

N-l lX N-l 



EC- 1 ) 



e=o 



N-l 



(i) n (x y y = (N-l) n J2(-lY 



N -n-l 



t-n 



{-lY{N-l) n {xy) n {l-xy) 



N-n-l 



Hence 

E[(S ab -b) n l {Sab > b} ] 

= (-l) n (N -l) n KN 2 



n-\-b— a— 1 / 



(l-xr-^-'il-xy) 



JV-n-l 



dx dy. 



Performing the change of variables (it, v) = (x, xy) in the foregoing integral immediately 
yields ( |8.30[ ). Formula ( |8.31[ l can be deduced from ( |8.26| l exactly in the same way. □ 

9 Link with high-order finite-difference equations 

Set A+/(j) = f(i + 1) - f{i) and A~/(i) = f(i) - f(i - 1) for any i <= Z and (A+)j = 
A+ o ■ ■ ■ o A+ and (A~) J ' = A" o ■ ■ ■ o A~ for any j e N*. Set also (A+) / = (A") / = /. 



j times j times 

The quantities (A+) J and (A - )- 7 are the iterated forward and backward finite-difference 
operators given by 



(A+)V(»)=E(-1) 



fe=0 



k 



fit + k), (A-ym = J2(-i) k KJ /(i - fc). 
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Conversely, f(i + k) and f(i — k) can be expressed by means of (A+) J f(i), (A y f{i), 
< j < k, according as 

f(i + k) = J2 ( k ) (A+ym, f(i -k) = E(-i) j ( k ) (a-)vw. (9.D 

j=0 j=0 

We have the following expression for any functional of the pseudo-random variable S a b. 
Theorem 9.1. We have, for any function f defined on E, that 



N-l N-l 

nns ab )\ = E + E ( A+ w) 

3=0 3=0 



(9.2) 



with 



N-l 

3 



u ] - a - l {\ - u) J "- 3 - 1 -v b - 1 (l - v) ly -' dudv 



\N-j-lb-l 



N-l 



v- 



It b ,3 = i-^YKN 2 
Proof. By (|9.1[), we see that 



N - 1 
3 



■D+ 



- u) N - 1 V^ +b - 1 (l - vf-3- 1 dudv. 



a—S a 



E[f(S ab )]=E\ l {Sab < a} E (- 1 

3=0 
Sab-b 

+ E\l {Sab > b} E 



3 



3=0 



N-l 

£(-ir'E 

3=0 
N-l 



{Sa. b <a} 



Sab - b 
3 

a — S a b 

3 



E E 

3=0 



{S ab >b} 



Sab - b 

3 



(A- )■>/(«) 

(* + y.f(b) 

(A-)V(o) 
(A + )V(&) 



which immediately yields ( |9.2| l thanks to ( |8.32[ l and l |8.33| l. 
Corollary 9.2. The generating function ofS a b is given by 

N-l N-l 

E(c s^ = c i- b .(i - i/cy + c b E 7 a\,-(c - !) j 

J=0 3=0 



(9.3) 



□ 



(9.4) 



Proof. Let us apply Theorem 9.1 to the function f(i) = ( l for which we plainly have 
(A+)V(i) = C(C - 1) J and (A~ )•*/(«) = C(l - 1/C) j - This immediately yields ( [9~4^ . □ 

Of special interest is the case when the starting point of the pseudo-random walk is 
some point x e Z. By translating a, 6 into a — x,b — x and the function / into the shifted 
function /(. + x), we have that 

N-l N-l 

E x [f(S ab )} = E[f(x + S a - X ,b- X )} = E K ,./. , - (A-)V(a) + E tf-«,6-xj 

J=0 J=0 



More precisely, we have the following result. 
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Theorem 9.3. We have, for any function f defined on E, that 



N-l 



JV-1 



E 



: [f(S ab )} = ]T P~ bJ (x) (A-yf(a) + KbjW (A + ) J /(fo) (9-5) 



3=0 



3=0 



where P ab ■ and P+ b < j < N — 1, are polynomials of degree not greater than 2N — 1 
characterized, for any fee {0, . . . , N — 1}, by 

(A-) k p- bj (a) = (A+) fe P+ d (b) = 6 jk , (A-) k P^ 3 (a) = (A+) k P~ bj (b) = 0. (9.6) 

Proof. By setting P^(a) = I^_ Xth _ x>j and P+ b / x ) = ^-xfi-^j' @ immediately 
yieldsl9.5l. By observing that P~ fc Ax) — (— l) J P+ b .(a+6 — x), it is enough to work with, 
e.g., Pabj- Coming back to the proof of Theorem 
write that 



9.1 



and appealing to Theorem 



8.6 



we 



pi.,(*) = E a 



{Sai,>b} 



Sab - b 



J 



N-l 



= E [ m WASab = b + m} 



N-l 



(N-l) 



V(-D T 



m=j 



m\ fN-l\ .fm + b-a + N -1 



J J \ m 



N 



K m {x) (9.7) 



where, for any m G {0, . . . , N — 1}, 



(N-l \ N-l 

Y\ [(x - a + k)(b - x + k)} ] /{b - x + m) = JJ (x - a + k) }~| (b-x + k) 
k=0 / fc=0 o<fc<N-i 



The expression K m (x) defines a polynomial of the variable a; of degree 2N — 1, so P+ b . 
is a polynomial of degree not greater that 2N — 1. 

It is obvious that K m {a - £) = for £, m e {0, . . . , N - 1}. Then, P+ M (a - I) = 
which implies that (A~) fc P+ h j(a) = for any k 6 {0, . . . , N — 1}. Now, let us evaluate 
P+ .(& + 1) for £ e {0, . . . , N - 1}. We plainly have that K m (b + £) = 0for£^ m and that 

+ b- a + N -l\ JN-1 



K e (b + £) = (-1) £ (N -l)\Nl 
By putting this into l |9.7[ l, we get that 

K b Ab + £) = 



N 



Next, we obtain, for any k € {0, . . . , N — 1}, that 

k 



1=0 



( A+ ) fcp ^( & ) = T,(- 1 ) k+ \ e ) p+ abA b + i ) = E(-!) fc+f 



k 

eJ\J 



(-i) 



j+k 



k-j 



k-j 
I 



= 5 



ijk- 



The proof of Theorem |9.3| is finished. 

Example 9.4. In the case where N = 2, ( |9.5| ) writes as 

E x [f(Sab)\ = Pa b , (x)f(a) + p- b>1 (x)A-f(a) + P+ {x)f(b) + P+ M (x)A+/(6) 



□ 
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with 



(x -b)(x-b-l)(2x-3a + b + 2) „_ , , (x - a)(x - b)(x - b - 1) 

P ab,l( X ) ~ 



(b-a)(b-a + l)(b-a + 2) 
(x - a) (a; - a + l)(2x + a - 36 - 2) 



P+ M W 



(b-a + l)(b-a + 2) ' 
(a; — a) (a; — a + l)(a; — b) 



(b-a)(b-a+l)(b-a + 2) ' * a6 - lw (6 - a + - a + 2) 

Below, we state a strong pseudo-Markov property related to time a ab . 
Theorem 9.5. We have, for any function f defined on Z and any n € N, that 

N-1 N-1 

E x [f(S aab+n )} = PabjW (A-)% Q [/(5„)] + ^ P+» (A+)% b [/(S n )]. (9.8) 

j=0 3=0 



In 



9.8), the operators (A ) J and (A+) J act on the variables a and b. 



Proof. Formula ( 9.8 1 can be proved exactly in the same way as (6.6 1: by setting g(x) = 
E x [f(S n )], we have that E x [f(S„ ab+n )] = E x [g(S ab )] . This proves <|9.8|l thanks to ^E\. □ 



Example 9.6. Below, we display the form of ( |9.8| l for the particular values 1, 2 of N. 
• For N = 1, l |6.6| l reads as 

E x [f(S aab+n )] = b ^E a [f(S n )} + ^E b [f(S n )} 



b — a 



which is of course well-known! This is the strong Markov property for the ordinary 
random walk. 

• For N — 2, l|6.6|l reads as 



= " 7,'"!^ ".t^" 



(b-a)(b-a + l)(b-a + 2) 

(b — a + l)(o — a + 2) 
(x - a)(x - a + l)(2x + a - 3b - 2) 

(6-a)(6-a + l)(fe-a + 2) 
(.x — a) (a: — a + 1) [x — b) 



E b [f{S n )\ 



(b- o + l)(6- a + 2) 
(a; - a + l)(x - &)(z - & - 1) 
(&-a)(6-a + l) 
[x — a)(x — &)(x — b — 1) 

(6-a + l)(&-a + 2) 
(x — a) {x — a + 1) (x — b — 1) 

(b-a)(b-a + l) 
(x — a){x — a + l)(x — b) 
+ (b-a + l){b-a + 2) 



A+E b [f(S n )} 

E a [f(S n )} 
E a -i[/(S„)] 
E b [.f(S n )} 



E b+1 [f(S n )} 



Now, we consider the discrete Laplacian A = A + o A = A o A + . It is explicitly 
defined by A/(i) = f(i + 1) — 2f(i) + f(i — 1). Let us introduce the iterated Laplacian 
a n _ (A+) JV o (A")^ = {Ar) N o (A+) 7V . We compute A N f{i) for any function / and any 
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i e Z: 



A\/(/) = (A-)-( >J-1I' + A (": )/(-+./)|(/) 

0<j,k<N \J / \ / 



AT /(N-t)/ 

HfEHH E 



=(-f)V0 



!/(< + *). 



By using the elementary identity Efc=(f- g )vo GDG-fe) = ( P | 9 )' we 9 et that 

k={-£)V0 x 7 v 7 fc=(-£)V0 



N-k-lJ U + N 



As a result, we obtain the expression of A N f(i) announced in the introduction, namely 



t=-N 



Example 9.7. Fix a nonnegative integer j and put fj(i) = for any i 6 Z, It is 
plain that, if k < j, (A+) fe /,(i) = (j) k £_*(*), (A") fe /iW - (j) fc /i-fc(< " 1) and if k > j, 
(A+) k f 3 (i) = (A-) k f,(i) = 0. Therefore, if 2k < j, A k f j {i) = (j) k (j ~ k) k f^ 2k (i - 1) and 
if 2k > j, A k fj(i) = 0. By using a linear algebra argument, we deduce that A N P = for 
any polynomial P of degree not greater that 2N — 1. As a byproduct, 

A N P+ b!j =A N p- b>j =0. (9.9) 

Now, the main link between time a a b and finite-difference equations is the following 
one. 

Theorem 9.8. Let <p be a function defined on £. The function $ defined on Z by <&(x) = 
E x [ip(S a b)] is the solution to the discrete Lauricella's problem 

A N $(x) ^OforxeZ, 

(A-) fc $(a) = (A-) V(a) for j £ {0, . . . , N - 1}, (9.10) 
t (A+) fc $(&) = (A+) k (p(b) for j 6 {0, . . . , N - 1}. 

Proof. By ( |9.5[ ) we write that 

N-l AT-1 



With this representation at hand, identities ( |9.6[ l and ( |9.9| l immediately yield equa- 
tions |9~TQ| |. □ 
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10 Joint pseudo-distribution of [r ab , X ab ) 

As in Section[7] we choose for the family ((X^) t > )e>o the pseudo-processes defined, 
for any e > 0, by 

X\ = eSy t / e 2N^ , t > 0, 

and for the pseudo-process (X t )t>o the pseudo-Brownian motion. In Definition [77T| we 
choose for / the interval (a, b); then r| = T^ b , Xj = X e ah and rj = r ab , X 1 = X ab . Set 
a e = [a/e\ and b e — \b/e~\ where [ - J and [ .] respectively stand for the usual floor and 
ceiling functions. We have T e ah = e 2N a a ^ be and X e ah = eS a ^b F _- 

Theorem 10.1. The following convergence holds: 

( T ab> X lb) (Tab,X ab ) 



where, for any A > and any /i e R, 

E{e- x ^ +i ^l {Tab<+oo} ) =e^«53^ 



AT t.- M s iV 
D fe( A '/ i ) , e i M 6 



(A) 



v d^(a 1 m) 

^ D(A) • 



k=l v ' fc=l 

In the foregoing formula, D(A), D,7(A,^), D^(A,^i) are the respective determinants 



1 <£x 



W-l p -¥>l 2 {/A/c(b-a) p - Vl WA/c(fe-a) 



H>2N 



e -^^fc(b-a) yN-1 



e -P2N 2 ^f>7~c(b~a) 



1 ipx ...if 



N-l p - Vl 2 &\/c(b-a) p - Vl 2 &\/c(b-a) 



1 Wk-l 

1 S 

1 ^2^ 



N-l --ipuS&X/cib-a) p - Vfc _i 2 ^A/c(b-a) 



¥>fe- 
riV-l 



<Pl 



Vk-1 



e - Vk + 1 ^X/c(b-a) e - Vk + 1 ^X/c(b-a) 



l/fa 1 e -V2 N 2 yiTc(b-a) e - V2N 2 -^xTc(,b-a) ^ 



- Vl 2 y\/c(b-a) N-l 



e - Vk ^X/c(b-a) ^N-l 



- Vk+1 2 y\/c(b-a) N-l 



p -V2iV 2 -(fA/ C (h-a) N-l 



1 ipi 

1 V'fe-l 



1 ¥k+l 



tpN-1 e - Vl 2 ^\Tc(b-a) 



<Pk-l 



^f+i 1 e-^^ l2 v^( b - Q ) e-^+ l2 v^( b - a )^. +1 



<^ e -^" 2 v^( b - a ) e-^" 2 v^( b - a )^ 2Ar 



1 <P2N 

In the two last determinants, we have put 5 



l\/c 



e - Vl 2 ^X/c(b-a) 



e -^ fc -l 2 {fA/c(fc-a) yN-1 
§ N-1 

- Vk+1 2 fITc(b-a) N-l 



e - V 2N 2 $ \/c(b-a) ^N-l 



In Theorem 10.1 we obtain the joint pseudo-distribution of (T ab ,X ab ) characterized 
by its Laplace-Fourier transform. This is a new result for pseudo-Brownian motion that 
we shall develop in a forthcoming paper. 
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Proof. By Definition |7 . 1 1 and Theorem [872] we have that 

i 



Urn E(e- x <» +i » x «* 1 {t *,< + oo}) 



{r ab <+oo} J = ™ + ^ - r ao J1{7 ^ < 

= lim IeT e _Ae2JVtT -' b - +i ' x£S '— b - 11 



{^a £ ,6 e <+00} 



2JY 



k=l 

2JV 



^E^^^^^i^i^)"" (10 - 1) 



Recall that L fe (z, £) = £>fc(z, ()/D(z) and that the quantities Z? and D k are expressed by 
means of the determinant 



W(ui, . . .,U 2 n) 

1 Ui u\ 



JV-1 „,&-a+JV-l „,fc-a+JV h-a+JV+1 



Ui U 



1 W2JV U 



2 

2iV 



?/ JV ~ 1 „fc-«»+^-l n .b-a+N b-a+N+1 

7\r ""J AT ">rj M tin 



,b-a+2JV-2 



,fe-a+2JV-2 



X 2JV "2AJ "2JV "2JV ■ • ■ "2JV 

By replacing the columns labeled as Cj, 1 < j < 2N, by the linear combinations Cj = 

£*=o(-l) J ' +fc if 1 < 3 < N, and C$ = Ei=jv(-l) J+fc (£)Cfc * ^ + 1 < j < 2N, the 

foregoing determinant remains invariant and can be rewritten as 



1 Ml - 1 (Ul - l) 5 



(ui - 1) 



JV-1 



1 U 2 N ~ 1 (U2N ~ l) 2 ■ ■ • («2AT - 1) 



AT-1 



b-a+N-1 n ,b-a+N-l 



h-a+JV-1 h-a+N-1 
a 2N a 



(tt-i) ui-^-^u-i) 



fc-a+2JV-2 
*1 



(«1 " 1) 



JV-1 



l 2N 



(u-1) u 2Ar 



b-a+N-1 



{u-lf ... M ^+ 2W - 2 ( U1 -1) 



\JV-1 



Then, by replacing the uj's by Uj{\, e), b — a by b e — a e and by using the asymptotics 

k 



(uj(X, e) - l) k ~ (-(pftfX/c) e k and Uj(A,e) 



b e -a e +JV-l 



e-s-0 

from l|4.5[l, we get that 



(c- Ae2 ") = w(ui(A, e), . . . , u 2N {\, e)) 



>0+ 



AT-1 



e->-0+ 



-ipi 2 {7A/c(6-a) 



coming 



1 A/V ) .: 



(-^jv 2 ^)^ 'e"- 1 e"*-^v v ' "' < 



c]e -^<7^/c(b-a) £ 



.JV-1 



e -¥>i 2 {/A7^(6-a) £ JV-1 



(-^2JV 2 ^A7^)e-^ 2 " 2 v y V^( fc - a ) e ... (-^ 2Af 2 ^A^ 

(A/c) (JV-l)/2 e JV(JV-l) D(A) ^ 



JV-1 



e -^2Jv 2 {fA7^(6-a) £ JV-1 



(10.2) 
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Similarly, by using the elementary asymptotics c 1/je — 1 ~ ins and e 1M£ ( fce a ?+ N i) — > 
e i/i(6-a) ; we obtain that 

e^ a D k (e~ Xs2N , e* 6 ) = e""W(ui(A, e), . . . , u fc _i(A, e), e* 6 , Ufc+1 (A, e), . . . 3 u 22V (A, e)) 

(10.3) 



£^0 + 



where Dfc(A,/z) denotes the determinant 

1 (pi . 



iV-1 _-«f.*_iWA/c(6-a) --^-^A/c^-a) 



1 (pk-1... <p k Z 
1 ^fe+i . . . <p^ + { v 



- Vk _^^/c(b-a) N-l 



p 2N ... e-^"^^-)^ 1 



1 ^ 2JV ... ^ e-^VVc(b-a) e -^^A/c(b-a) 

By putting ( fTQT2| l and ( p~073] > into ( fTQTTj l, we derive that 

F ( -^ ab +hiX ab \ _ D fc( A ^) 

-H{r ab <+oo} J - 2^ D(A) ' 

It is plain that (A, fi) — e llJ,a T)7 (A, /i) + e'^D^ (A, /x) which finishes the proof of Theo- 
rem corn □ 

Theorem 10.2. The following convergence holds: 

Xth — > X„h 



where, for any (jgl, 



E 



with 



AT-1 JV-1 
J=0 j=0 



b \ a J x fk + N — l\f —a 



abj U - a 7 7'! 



E 



?! \ A; 



6 — a 



a6j V b - a 7 ?! 



E 



-a r (-^ "^VHJV-A/ b 



J fe=0 



b- 



Moreover, 



¥{r ab < +00} = 1. 
Proof. By definition [77T] and by ([9T4[>, 

E(e i ^-l {rab<+oo} 
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(N-l N-l 
^ E - e ^ £ y + ei " £be E - i y 
3=0 3=0 

N-l N-l 

e^ a lim V/" 6 ..(l-e-^'+e^ 6 lim V 1+ b Ae^ s - l) j , (10.4) 
Concerning, e.g., the quantity 7+ be we have that 



u-^-^l - u) N - 1 v s+b '- 1 (l - v)*-!- 1 dudv. (10.5) 

v+ 

By performing the change of variables (u, v) = (x, xy) in the above integral and by 
expanding (1 — xy) N ~^ x as 

(1 - xy)^- 1 = [(1 -x) + x(l - y)}^- 1 

N-j-l , . _ x 

= E ( * Jci-^^-^a-i/)^'-*- 1 , 



fc=0 

we get that 



u- a °- x {\ - ^^"V'+^-^l - v)^- 1 dud?; 
p 

i /.i 

j+bs-as-l^ _ ^N-lyj+be-l^ _ xy jN-j-l dxdy 

oJo 

N -l- 1 /AT A 1\ /•! /-I 



E f ^ /'/-^^-"(l-^^^Mx / 1 ^- 1 (i- ?/ )^- fc - 1 dj/ 

l._n V k J Jo Jo 

( N - j - 1\ (b e - a e + N - k - 2)\{k + N - 1)! (j + b £ - l)\(N -j-k-l)\ 



k J (b e - a E + 2N - 2)! (6 e + JV-fc-l)! 



fc=0 

N-j-l 



(N-l)l(N-j-l)l(j + b s -l)\ ^ (k + N-l\ (b e -a e +N-k-2)\ 



(b e - a e + 2N - 2)! ^ V fc / (6 e + iV - - 1)! 

By putting the asymptotics 

(j + b e -l)\ 1 (e\ N -3-k 



(b e + N-k-l)\ (b s +N -k-l) N _j_ k £ ^o+ V& 
(b E -a e + N-k-2)! _ 1 / E 



(b e - a £ + 2N - 2)! (6 £ -a e + 2A^-2) A , +fe £ ^o+ V&-a 
into l |10.6| l, we obtain that 

u- a '-\l - u) N - 1 v 1+b '- 1 (l - v^- 3 - 1 dudv 

p 

(N -l)l(N - j -1)1 jj^-VT 1 fk + N-l\f b 



e 



E 
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Next, using the asymptotics 

JV-l 



l[[(k-a e )(k + b e )}^ o+ (-1) 

k=0 



N 



(ab) N 

r 2N : 



expression ( |10.5[ l admits the following asymptotics: 

{-iy+ N v ( a \ NN ^ 

k=0 



Then, we see that the second limit lying in ( |10.4[ ) tends to 



k + N -1 
k 



b — a 



b — a 



E 

3=0 



J! 



N-j-l 

E 

L fc=0 



k + N-1 
k 



b — a 



In the same way it may be seen that the first term of the sum lying in l |10.4[ l tends to 



b — a 



N JV-l ,. 



E 

3=0 



(i^op 



'N-j-l 

E 

, fe=0 



k + N-1 

k 



b- 



Finally let us have a look on the pseudo-probability P{r ab < +00} . We have that 



L ab,0 



, \N JV-l 

b-a J E 
b 



k=0 
A N-l 



k + N-1 
k 



b — a 



(b-a) 



2JV- 



rE 

k=0 



k + N-1 

k 



(-a) k (b-a) 



N-l-k 



f,A r 



(b-a) 2N - 
1 



T E 



0<fc<iV-l 
0<£<iV-l-fc 



k + N-l\fN-l-k 



k 



(-a) k+t b 



k+e tjv-i-fc- 



(b-a) 



2JV- 



JV-1 

rE 



E 

fc=0 



k + N - l\fN - 1 - k 
m — k 



(-a) m b 



rn l2AT— 1—771 



By using the elementary identity J2k=o ^V^t-k") = ^ l+P f +1 ) which comes from 
the equality (1 + x)~ p (l + x)~ q = (1 + x)~ p ~ q together with the expansion, e.g., for p, 

(1 + x)-p = E^ot-^CTV' we 9 et that 



E 



k + N-l\fN-l-k\ (2N-1 



As a byproduct, 



Similarly, 



JV-l 



L ab.O 



(b-a) 



1 y 



m— 



2 JV-l 



(b-a) 



1 y 



m=N 



and we deduce that 

P{r Qb < +00} = I" w + 1+ 



- k 



2N - 1 
m 

2N - 1 
m 

2 JV-l 



(-a) m b 



mi^N—l—m 



(b-a) 

The proof of Theorem 1 10.2 1 is finished. 



1 y 

A2N-1 



(-a) m b 



(2N-1 
\ m 



2JV-1-T 



(-a) m b 



2JV-1-T 



□ 
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Corollary 10.3. The pseudo-density of X ab is given by 

JV-l N-l 

¥{x ab e dz}/dz = £ sP(z) + J2 i- 1 )^ ^\z). 

In particular, 

P{r.-<7-+} = I- 6i0 , P{r+<r-}=I+ )0 . 

This result has been announced in H13B without any proof. We shall develop it in a 
forthcoming paper. 
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A Appendix 

A. 1 Lacunary Vandermonde systems 

Let us introduce the "lacunary" Vandermonde determinant (of type (p + r) x (p + r))-. 



W pqr (ui, . . . ,u p+r ) 



1 Ui 



p— 1 P+q 



p— l „,p+q 



p+q+r — 1 



p+q+r — 1 



1 Wp+r • • • "p+r "p+r ■ ■ ■ "p+r 

We put s := sq{u\, . . . , Up +r ) = 1 and, for I e {1, . . . ,p + r}, 

se := se(ui, . . . ,u p+r ) = ^2 u il ---u it . 

l<i±<---<il<p+r 

We say "lacunary" because it comes from a genuine Vandermonde determinant where 
the powers from p to (p+q— 1) are missing. More precisely the determinant W pqr (ui, . . . , 
u p+r ) is extracted from the classical Vandermonde determinant V p+q+r (ui, . . . , u p+q+r ) 
of type (p + q + r) x (p + q + r) by removing the q last rows and the (p + l)st, (p + 2)nd, 
. . . , (p + q)th columns. We decompose V p+q+r (ui, . . . , u p+q+r ) into blocks as follows: 



1 



Ml 



1 Wp_|_ r . . 
1 Up+ r +l 

1 ^ p + q -\-r • ■ 



,P-1 



p-1 
^p+r 



p+q-1 



,p+q 



P+q-1 P+q 
"p+r a p+r 



.p+q+i — 1 



p+q+l — 1 



P+r 



p—1 p 
U p+r+l U p+r+l 



u p+q+r u p+q+r 



p+q-1 p+q 
"p+r+1 "p+r+1 



p+q-1 p+q 
U p+q+r u p+q+r 



p+q+r-1 
i p+r+l 



p+q+r-1 
l p+q+r 



By moving the r last columns before the q previous ones, this determinant can be rewrit- 
ten as 



1 

1 l^p-yr 
1 U p+r+ l 



p-1 



p-1 

*p+r 



,p+q 



p+q 

l p+r 



v p-i v v+q 

a p+r+l a p+r+l 



1 U p+q+r ... U 



P-1 

p+q+r "p+g+r 



p+q 



U 



p+q+r-1 



U 



p+q+r-1 
p+r 

p+q+r—1 
p+r+1 



U 



p+q+r—1 
p+q+r 



l p+r 



U 



p+r+1 



U 



p+q+r 



u; 



p+q-i 



u: 



p+q-i 

p+r 

p+q-1 
1 p+r+1 



U 



p+q-i 

p+q+r 



By appealing to an expansion by blocks of a determinant, it may be seen that W pqr (ui, 
...,u p+r ) is the cofactor of the "south-east" block of the above determinant. Since 
the product of, e.g., the diagonal terms of this last block is Up + r-+i u p+r+2 ' ' ' u p+ 9 q +l> 
the determinant W pqr (ui, . . . , u p+r ) is also the coefficient of «p + r+i M p+r+2 ' ' ' u p+g+r i n 
V P +q+r(ui, ■ ■ ■ ,u p+q+r ). Now, let us expand V p+q+r (ui, . . . ,u p+q+r ): 

Vp+g+rim, U p+q+r ) = Yl ( U 3 - U i) = III Il 2 IT3 

l<i<j<p+q-\-r 
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with 



iii = n k - ^ 

l<i<j<p-\-r 



n 2 
n, 



n {uj-ui), 



l<i<p+r 
p+r+l< j < p+q+r 



n (uj-ut). 

p+r+l<i<j<p+q+r 



We have that 



o+q+r 



p+q+r p+r 



n = n [(u j -u 1 )---(u j - Up+r )}= n 



j=p+r+l 



j=p+r+l k=0 



E(_l)g(p+r)-ki k q , ? , fc i ■■■ V kq 

V 1 ) Op+r-k! *p+r-k q « p + r +i « p + g + r 



0<fci,...,fc,<p+r 



and 



U 3 — Vq( u p+r+l, ■ ■ ■ ,U p+q+r ) — ^ t{s) u y+r+l ' 



2 e(4...,^X +r+1 ' 



■ W 



0<<i,...,/ 9 <9-l 
£i ,. .. ,e a all different 



The symbol & q in the above sum denotes the set of the permutations of the numbers 
1,2,..., q, e(c) is the signature of the permutation <r and e(£i . . . ,£ q ) e {— 1, +1} is the 
signature of the permutation mapping l,...,q into £\...,£ q . The product \\ 2 ]J 3 is given 
by 



( _1)<Kp+i) (-lp + - +k ^(l 1 ...,l q )s p+r ^ kl --- 3p+r ^ 



<, a p+r+l a p+q+r- 



0<k 1 ,...,k q <p+r 

o<e 1 ,...,e q <g-i 

t-t ,...,C a all different 



For obtaining the coefficient of u p+r+1 u p + r+2 ■ ■ ■ Up+g+r, we on ly keep in the foregoing 
sum the indices ki , . . . , k q , t\ , . . . , £ q such that k\ + £ i = p, ki + £2 = p + 1, . . . , k q + £ q = 
p + q - 1 and < k\ , . . . , k q < p + r, < l\ , . . . , l q < q — 1, the indices £\,...,£ q being all 
distinct. This gives that 



n 2 n : , 



c{£\ £q)s£ 1+r S£ 2+r - 1 ■ ■ ■ S£ q+r - q+ l. 



fori£{l, ...,<,}: 

(i-i l jVO^ti <(p + i-l)A(g-l) 

e 1; ...,e q all different 



Finally we can observe that the foregoing sum is nothing but the expansion of the 
determinant 

S T S r —\ . . . S r — q -\-\ 

°r+l &r • ■ • S r — q +2 

^pqr 

&r+q— 1 &r+q— 2 

As a result, we obtain the result below. 

Proposition A.l. The determinant W pqr {u\, . . . , u p+r ) admits the following expression-. 

W pqr {u 1 ,...,u p+r )= Y[ ( U j~ U i) 

l<i<j<p+r 
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s r-q+l{ u l, ■ ■ ■ j u p+r) 
S r -q+2(ui, • • • , U p+r ) 

S r {ui , . . . , tip_|_ r ) 



Let now W pqr ^ I ' ' ' ' "' p+r I be the determinant deduced from W pqr (ui, . . . , u p+r ) 

\OLl , . . . , Qfp+r / 

by replacing one column by a general column ai, . . . , a p+r , that is, the determinant 
given, if < £ < p — 1, by 



s r («i, . . • ,Up+ r ) s r _i(ui, . . . ,u p+r ) 
s r+ i(ui, . . . ,u p+r ) s r (ui, . . . ,u p+r ) 

&r+q— 1 

(til, • • ■ ,Up+r) S r+q - 2 {ui, . . .,Up +r ) 

U\ , . . . , tip-|_ r 



1 til 



1 ti 



«1 ti x 



£-1 ,/+! 
ip_|_ r Ctp+r u p+r 



and, ifp + £jr<£<p + <7 + r — 1, by 



1 tii 

1 U p -^- r 

We have that 



u p ~ l u p+q 



p-i p+? 

"p+r "p+r 



p— 1 P+g 
ti, Ml 



p-i „,p+g 



ai tif 1 



,p+g+r— 1 



,p+g+r— 1 



x p+r 



p+g+r-1 



u 1 - 1 a . u e+1 

"p+r "p+r "p+r 



u: 



p+g+r-1 
p+r 



A ' ' ' '" P+r ) = y"^afcWpgr.fcl(tti, . . . ,Ufe_i,U fe+ i, . . . ,tip+ r ) 

yai, . . . , a p + r y j 



where W pqr ^i{ui, . . . , Uk-i,Uk+i, • • ■ , ttp+ r ) is the determinant given, if < £ < p — 1, by 



1 ui 



4- 1 o 



. . . u' 



p-i „,p+g 



ti 



p+g+r-1 



Ufc-l • 







< + -\ ■ 


p-1 

• <-l 


u p+q 
"fe-i ■ • 


p+g+r- 
• "fc-1 


1 


Wfc • 




1 




p-1 

• Ufc 


"fc 


p+g+r- 

• "fc 


1 


ttfc+1 • 









p-1 

• "fc+1 


"fc+1 ' ' 


p+g+r- 
• "fc+1 


1 



.£-1 



""p+r ti p+J , 



1 tip_j_ r • 

and, ifp + <7<^<p + <7 + r — 1, by 



,t+i 



1 tii 



1 tifc_i ... 
1 tifc 



ti?" 1 u p+q 



p—1 p-\-q 
U p+r "p+r 



4- 1 o u[ +1 



p+g+r— 1 
p+r 



11 



p+g+r-1 





u p+q 

"fc-1 ' 







4t\ .. 


p+g+r 
• "fc-1 


r 1 


"fc 




1 




p+g+r 

• "fc 


p-i 
-fc+i 


"fc+1 ' 







t/ +1 

"fc+1 • • 


p+g+r 

• "fc+1 


p-i 

'p+r 


1l P+q 
u p+r 


U /_1 
• • u p+r 





u e+1 


p+g+r 
u p+r 



1 tifc+1 . . 

1 u p+r ... u 

In fact, W pq r t ki{ui, . . . , tifc-i, tifc+i, • • • , u p + r ) is the coefficient of wf i n W pqr (ui, . . . , ti p + r ). 
Let us introduce s fej0 := Sfc,o(wi, . . . , Ufe-i, Ufc+i, • • • , «p+r) = 1, «fc,m := Sfc,m(wi, . . . ,tifc-i, 
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Ufe+i, . . . ,u p+r ) = for any integer m such that m < — 1 or rn > p + r and, for m £ 
{l,...,p + r-l}. 



Sk,m '■— Sk,m(ui, ■ ■ ■ , Ufc_l, Ufc+1, . . . , U p+r ) — ^ 



l<il< ■ <i m <p+T 
il,...,im^* 



We need to isolate itfc in rii<;<j<p+r( M j — u i) anc ^ ^pgr- First, we write that 



l<i< j<p+r 



l<i<3<p+r 



l<i<p+r 



P+1 1 



= (-l) fe ~ 1 J] (uj-Ui)x J] (-l) m ^_ m . 1 <. (A.1) 



l<i<j<p+r 



Second, by isolating itfe in s m according as s m — s k ^ m + "fe Sfe m-i< we !J e t that 



Spqr 



Sfe,r + U fe Sfe, r -i 



Sfe,r-1 + w fc s fc,r-2 



Sfe,r+g-l + Ufe S k .r+q-2 Sk.r+q-2 + Ufc Sfc. r + g _3 



Sk,r— g+1 *T 'Wfe Sk.r—q 
s k,r-q+2 + u fe s fc,r-ij+l 



By introducing vectors V^n with coordinates (sfe, m , Sfc, m +i, . . . , Sk, m +q-i) (written as a 
column), this determinant can be rewritten as 

Spqr = det(V r + U k V r -l,V r -\ + U k V r -2, • • • , VV-g+1 + U k V r - q ). 

By appealing to multilinearity it is easy to see that 



Spqr — det(t? r , V r -i, . . . , V r -q-\- n -\-l, V r - q -\-n-l, . . . , V r - 



q) a k- 



(A.2) 



n=0 



Now, let us multiply the sum lying in < |A. 1 1 > by l |A.2| l: 

/p+r-l W q _ 

(-i) m s k ,p+r- m -iu% J I ^det(V r ,Vr-l,...,V r - q+n+1 ,V r - q+n -l,...,V r -q)ul 
\ m=0 ) \ n=0 , 

p+q+r-1 / (p+r-l)M 

= Z2 [ X! (- l ) ms k,p+r-rn-l 
e=o \m=0V(e-q) 



det{V r: V r -1, ■ ■ ■ , Ve+r-q-m+l, Ve+r-q-m-1, ■ ■ ■ , Vr-q) u l- 



(A.3) 



Recalling the convention that s k _ m = if m < — 1 orm >p + r, the coefficient of u e k in 
( |A.3p can be written as 



(— l) m Sfe,p+ r - m -l dct(V^, V^- 1) ■ • • , Vt+r-q-m+l> Ve+r-q-m-1, ■ ■ ■ j 

m— f — q 

= ( _ 1^ l) m+9s /£,p+q+r-£-m-l det(V^, . . . , V r -m+l, Vr-m-li ■ ■ ■ i Vr-q) ■ 

m=0 
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In this form, we see that the coefficient of u e k in I A.3 1 is nothing but the product of (— 1) £ 
by the expansion of the determinant 



S, 



pqr,kt 



s k.r+l 



°fc,r— 1 
Sk,r 



Sk,r+q-l 8k,r+q-2 
Sk,p+q+r~£~l Sk,p+q+r-l-2 



Sk.r~q 
Sk,r-q+l 

Sk,r — 1 
Sk,p+r-£-l 



Proposition A. 2. The determinant W pqr ,kt(ui, ■ • ■ , Ufe-i, Uk+i, ■ ■ ■ , u p + r ) admits the fol- 
lowing expression: 

W pqr ,kt{Ul, ...,Uk-l,Uk+l,...,V, p+ r) = (-l)^ -1 Yl ( u j- u i) 

l<i<j<p+r 
i,j^k 

Sk,r-q(ui, ■ ■ ■ , Ufc-1, ltfc+1, ■ ■ ■ , U p+r ) 



Sfc,r(ui, ■ • ■ , Ufc-1, Uk+1, ■ ■ ■ , Up+r) 
Sfc,r+l(«l 5 ■ • ■ ,Uk-l,Uk+l, ■ ■ ■ ,U p+r ) 



Sfe,r-5+l(«l) ■ ■ ■ i Wfc-1) "fc+li • • • j Up+r) 



Sfe.r+g-l^l, . . . , lifc-l, Ufe+1, ■ ■ • , U p+r ) . . . S k , r -l(ui, . . . , lifc-1, Ufc+lj ■ ■ ■ , U p+r ) 
s fc,p+g+r-£-l(wi, . . . , Ufc-1, Ufc+1, . . . , U p+r ) . . . Sfc,p+r-£-l («1 , ■ ■ ■ , Ufc_l,«k-fl, . . . , ttp+r) 

As a consequence of Propositions |A. 1 1 and |A. 2| we get the result below. Set 

Pk :=Pk(ui,...,u p+r ) = Yl (uk-Ui). 

l<i< P +T 

Proposition A.3. Let p,q,r be positive integers, let u\,... ,u p+r be distinct complex 
numbers and let a±, . . . , a p+r be complex numbers. Set I = {0, . . . ,p— 1}U {p + q, . . . ,p + 
q + r — 1}. The solution of the system Yltei u l Xe = ak ' 1 < k < p + r, or, more explicitly, 



X + UlXi+"-+Ui 1 Xp-i + U^ x p+q + ■ ■ ■ +u 



p+q , 



p+q+r—l 



'P+q+r— 1 — Oi\ 



. ■ , p— 1 , p+q , , p+q+i — 1 

k ^O^Wp-t-r X\+- • • • ~rUp_(_ r X p —\ + U p _^_ r X p J rq -\- ' ' ' '+U p j rr X p +q-\-^ — \ Q^p-j-r 



is given by 



I 1 "li+p+r— 1 P+ r c I \ 

<->pqr,k£\ U l> ■ ■ ■ ; "fc + 1; • ■ • ; Up+r) „ j ( K A\ 

%i — c, i \ / . Oik i \ , -c fc i . 

6 pqr (Ui, . .. , Up+r) ^ Pk(Ux,..., U p+r ) 



Proof. Cramer's formulae yield that 

W pqr J 



Ul , ■ ■ • j U p -\-r 
Oil , ... , CKp_(_r 



p+r 



- / a k 

W pqr {Ul, . . . ,U p+r ) J W pqr (Uii . . . ,U p+r ) 

By using the factorisations provided by Propositions |A. 1 1 and |A. 2 1 namely 

W pqr (ux, . . . ,Up +r ) 

= Vp+ r (ui, . . . , U p + r ) S pqr (ui 7 . . . , U p + r ) 

= (-l) P+r_fc Pfe(ui, ■ • ■ ,U p+r ) Vp +r _i(ui, . . . ,Uk-l,Uk+l, ■ ■ ■ , U p+r ) Spqr(ui, . . . ,U p+r ), 

W pqri ki(ui, . . . , v,k-i,Uk+i, ■ ■ ■ , u p+r ) 

= + _1 ^+r- l(ui, . . . ,Ufe_l,Ufe+l, . . . , U p+r ) S pqr ,kl(ui, . . . , Ufc_i, Ufe+i, . . . , u p+r ), 

we immediately get ( |A.4[ I. 



W pq r t k£{Ul, . . . , Uk-lfUk+l, ■ ■ ■ , U p+r ) 



□ 
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A. 2 A combinatoric identity 

Lemma A. 4. The following identity holds for any positive integers ct,fi,n: 



k (n\ (k + a)\ 



k=0 

It can be rewritten as 



kj(k + p)\ (0 + n)\ 



(P -a + n- 1)„. 



f{-D k ( n )/( k+p )= ^ a /p + 

^ ' \k)'\k + a) p-a + n'\ a 



if a > (3, 
if a < /3. 



Proof. Suppose first that a > j3. Noticing that 

(fc + oQ! _ d"-' 3 
(k + /3)l ~ dx a -' 3 

we immediately get that 



(x k+a ) 



x=l 



n\(k + a)\ d a -P 



S ( 1)k {k)(k + py. dw B 1)A 



We expand the last displayed derivative by using Leibniz rule: 



dm 



\x a {l-x) n ) 



x=l 



&x a -P 



(x a (l-x) n ) 



x=l 



E 

k=0 



a-0\ d 



j«— 8— k 



k / dx«-^- fe 



dx fc 



x=l 



Since (dVdx fc ) ((1 - a;)") | x=1 = {-l) n S nk n\, we have that (d a -P /da;"-' 9 ) (x a (l - x) n ) \ x=1 
= if a — /3 < n, and, if a — /3 > n, 



dx a ~P 



(x a (l-x) n ) 



fa - B\ d Q_/3_n 
= (-l)"n' - 



j;=l 



(-1)" 



n / da;"-"" " 
a!(a-/3)! 



x=l 



{l-xf- a - l dx, 



{a - - n)\(P + n)\ 

which coincides with the announced result. Second, suppose that a < f3. Noticing that 

(k + a)\ _ 

JkTpy-~(P 

we get that 

^kf n \ (k + a)\ _ 



_i f 1 T k+» 

a-iy.Jo 



£(-d* 



k=0 



kj{k + p)\ {p-a 
1 



(J3-a-l)\J 
which coincides with the announced result 



x a {l- X y 



k) xk j x<X { l ~ x f~ a ~ lAx 
-a+n-i dx = a\(P-a + n- 1)! 



(/?- a- !)!(/? + n)! 



□ 
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A. 3 Some matrices 

Let a,/3 eN such that a > (3 > N and set 

f 3 + 



A = 



N 



B = 



0<i,j<N-l 



P 

i + N 



0<i<N-l 



with the convention of settings Q) = if i > j. These matrices have been used for 
solving systems ( |8.13[ ) and ( |8.14[ ) with the choices a = b — a + N — 1 and /3 = N — a — 1. 
The aim of this section is to compute the product of the inverse of A by B, namely 



A 1 B. For this, we use Gauss elimination. The result is displayed in Theorem A.8 The 



calculations are guite technical, so we perform them progressively, the intermediate 
steps being stated in several lemmas (Lemmas |A.4| |A5] |A.6| and |A. 1) . 

Lemma A. 5. The matrix A can be decomposed into A = LU^ 1 where the matrices U 
and L are given by 

r j\ {j + o)n 



U = 



L = 



h 



{i<3} 



j + a 



i j (i + a)jv 



0<i,j<N—l 



{l ^ } \i + Nj (j + a-N)j 



0<i,j<N—l 

The regular matrices U and L are respectively upper and lower triangular. 

Proof. We begin by detailing the algorithm providing the matrix U. Call Cq°\c[°\ 
C < fi ) _ 1 the columns of A, that is, for j e {0, . . . , N - 1}, 



C 



(0) 



j + a 
i + N 



0<i<N-l 



Apply to them, except for Cq°\ the transformation defined, for j 6 {1, . . . ,N — 1}, by 



j + a 



C (1) = C (0) - 
3 3 j + a 



N C i-v 



i Cn-i are the columns of a new matrix L\. Actually, this transforma- 
tion corresponds to a matrix multiplication acting on A: L\ = AU\ with 



1 





Simple computations show that 



a+l 
' a+l-N 

1 







g+2 
" a+2-N 

1 



Li 



Sjo + %>i} 



j + a 



j + a 
j + a- Nj\i + N 

(*)jAl 









a+N-1 
a-1 
1 



0<i,j<N-l 



0<i,j<N-l 



K i + NJ (j + a-N) jA1 _ 

Next, apply the second transformation to the columns of L\ except for and C^~\ 
defined, for j e {2, . . . , N - 1}, by 

j + a 



C (2) = C (1) - 



j + a - N 

j + a q(0) 



(j + a)(.l + a-l) 



C 



(i) 



= c (o) _ 2 

j j + a-N ^i- 1 ^ (j + a- N)(j + a - N - 1) ^~ 2 ' 
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The Cq°\ C[ x \ C 2 2 \ . . . , cffi-i are the columns of a new matrix L 2 = AU 2 where 



U 2 



g+l 
' a+l-N 

1 









(a+2)(g+l) 
(a+2-N)(a+l-N) 

a+2 



a+2-N 



(a+N-l){a+N-2) 
(a-l)(a-2) 

_9 a+JV-l 
a — 1 



Straightforward algebra yields that 



j + a - TV 

i + a\ (») 3 A2 

* + NJ (j + a- N) jA2 



5 n + l {i > 2} 



i(i-l) 



(j + a - iV)(j + a - TV - l)/ V« + N 



j + a 



0<i,j<N-l 



0<i,j<N-l 



This method can be recursively extended : apply the fcth transformation (1 < k < N — 1) 
defined, for j e {k, . . . , N — 1}, by 



r (fe) _ r (fe-i) _ j + a r (k-i) 
* j + a-iV 



c (fe-2) _ 2 j + a c (k-2) + (j + a)(j + a-l) c ( k - 2 ) 



j + a- N 



(j + a- N)(j + a- N - 1) ^ 2 



1 y j + a - N 



+ ■■■ + (-1)* 



+ + r (o) 

2j (j + a- N)(j + a- N- 1) J '~ 2 

(j + o:)(j + a - 1) • • • (j + a - fc + 1) 

V (j + " - N )U + a - N - 1) ' ' ' (j + a - k - N + 1) 



r (o) 



/" fcA (j + a)jv r (0) 
£j(j + a-£) N i~ l 



V r-iv'-'f fc ^ (i + g)jz* 
4^ j Kj-eJti + a-N)^ 



r (o) _ / k \ U + a )N r (o) 



In particular, 



The C^ 0) , Cf \ Cf \ . . . , C^ J , C£> v Cfi^ are the columns of the fcth matrix L k = AU k 



•,(k) n (k) 



-.(*) 
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with Uk = [U' k : U'l ] where U' k is the matrix 



a+l 



(a+2)(a + l) 



a+l-N (a+2-N){a + l-N) 



9 Q+2 
^ a+2-N 



V L > \k-lJ (a+k-N-1) 

I l\k-2(k-l\ (a+fc-1) 

V L ) \k-2l (a+fe-JV-1) 



/- ~\\k-3(k-l\ (a+fc-1) 
V \k-3J (a+k-N-1) 



(" + !) 



(a+l-N) 
("+2) 



(a+2-N) 
("+3) 



(a+3-N) 



(V) 



(g+fc-l)(g+fc-2) 
(a+fc-Af-l)(a+fc-JV-2) 

_(k-l\ a+fc-1 
VI/ a+fc-JV-l 

(V) 





and XJ'l is the matrix 



7 lUyfc'i (a+fc)-(g+l) 
V ^ Vfc/ (a+fc-AT)---(a+l-Af) 



(2) (a 



(a+fc)(a+fc-l) 



(a+fc-AT)(a+fc-Ar-l) 
, a+fc 
a+fc-Af 







(-l) fe 



(a+fc+l)-(a+2) 
,fc/ (a+fc-AT+l)...(a+2-Af) 



fc/fc\ 



(2) 



(a+fc+l)(a+fc) 
(a+fc-AT+l)(a+fc-JV) 

a+fc+1 
VI/ a+fc-. 



+fc-N+l 
(0) 



/ i^m (a+Af-l)-(a+Af-fc) 
V L > \k) (a-l)...(a-fc) 



/fc\ (a+jV-l)(a+jV-2) 
\2) (a-l)(a-2) 

_ (k\ g+N-1 
VI/ a-1 



(0) 



The matrices UL and UL' can be simply written as 



U' k = 



j\ U + oQjv 

i) (i + a) N 

k \ ( j + <*)n 



0<i<N-l 
0<j<k-l 



j -ij (i + a) 



N 



0<i<N — l 
k<j<N-l 



Clever algebra yields that 



Lh = 



(J) 



jAk 



j + a 

i + Nj (j + a- N) jAk 



0<i,j<N-l 



(A.5) 



We shall not prove I A.5 1, we shall only check it below in the case k = N — 1. 
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We progressively arrive at the last transformation which corresponds to k = N — 1: 



c (N-i) = c (o) _ (N - l\ N + a-l n{0) ^(N-l^ (N + a-l)(N + a-2) n{0) 



a-1 



N-2 



(a-l)(a-2) 

r-1 ) N - X ( N ~ l \ ( N + a - 1 )( N + a- 2 )---( a + l ) MO) 
[ ' \N-lJ {a-l)(a-2)---(a- N+l) 



JV-3 



JV-1 



N-l\ (N + a-l)i (o) 



1 ^ ~ A ( N + Q - 1 )n-i-i Mo) 



The Cq°\ C} 1 , C2 , . . . , Cj^_ 1 1 ^ are the columns of the last matrix given by £at_i 
-Af^v-i where 



(-1 



i) (i + a) N 



0<i,j<N-l 



Formula ( A.5| > gives the following expression for L^-i that will be checked below: 



Ln- 



j + a\ 



i + NJ (j + a - AO j 



0<j,j<JV-l 



Hence, by putting L = L^-i and U = Un-i, we see that L is a lower triangular matrix, 
U is an upper triangular matrix and we have obtained that L = AU. 

Finally, we directly check the decomposition L = AU. The generic term of AU is 



Vr \v+ k ( k + a \( j \ + 

\i + Nj\kJ(k + a) 



N 



Observing that 



K i + NJ (k + a) N 
this term can be rewritten as 



k + a\(j + a) N _ (k + a- N\ (j + a\ fi + N 



N 



N 



(-iy 



n r \ n 



k=0 



'j\fk + a-N 



A.4 



Its 



The sum J2i=o(~^) k (i)( k+a i N ) can ^ e explicitly evaluated thanks to Lemma 
value is (— ). Therefore, we can easily get the generic term of L and the proof 
of Lemma IA. 51 is finished. □ 

Lemma A. 6. The inverse of the matrix L is given by 



U i+ S i {j + N)l(i + q-N) i+1 
(-1) J {i>j}^vr- ttt^ — \ 



0<i,j<N-l 



Proof. We simplify the entries of the product 



0<i,j<N-l 

(iV+j, (j + + a - N) i+1 

( ^ Hi>» jl(i _ mi + a)j+N+i 



(A.6) 



0<i,j<N-l 
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The generic term of this matrix is 

N-l 

E 



k=0 



k + a 



= 1 



i + NJ (k + a - N) k 

+ N)l(a - N)\ 
+ N)\(a - N - 1)! 



/ 1V -+*1 {j + N)l(k + a-N) k+1 
( L} ^> j\(k-j)\(k + a) j+N+1 



EC- 1 ) 



3 +k_ 



k=j 



(k + a-j-N-l)\ 
(i — k)\(k — j)\(k + a — i — N)\' 



The last sum can be computed as follows: clearly, it vanishes when i < j and it equals 1 
when i = j. If i > j, by using Lemma |A.4| 



y r lV 7+fc {k + a-j 



N-l)\ 



k=j 



(i - k)\(k - j)\(k + a - i - N)\ 



(i-j) 



EC- 1 ) 



(i-j)! 



k=j 
i—j 

EC- 1 ) 



j+fc ^ - A (fc + g-j-JV-l)! 

fc-jy (fc + a-i-iV)! 



1 1-3 

1 N - -^k l -J 



fe=0 



(fc + a - JV - !)<_,•_! = 0. 



As a consequence, the entries of the product lA.6l are Sij which proves that the second 
factor of l A. 6 1 coincides with L^ 1 . 



a 



Lemma A. 7. The matrix L 1 B is given by 



L~ B = 



+ a — — \\ fi + a 
N 



y NJ\ a- 0-1 
Proof. The generic term of L~ 1 B is 

{j+N)\(i + a-N) i+1 f 



0<i<N-l 



N-l 



Ec- 1 )^' 1 ^} 



3=0 



- + a) j+N+1 \]+N 



0\(i + a-N)\ 
(i + a)\(a - N - 1)\ 



EC- 1 ) 



i+j 



J=0 



(i - j + a - N - 1)! 
jl(i-j)KP-j-N)\ 



0!(i + a- N)\(i + a-0-l)\^ f 1 f i\f i — j + a — N — 1 



+ a)\(a - N - 1)! 



Ec- 1 )' 



iy \ i + ol — — l 



(A. 7) 



By performing the change of index j n- i — j and by using Lemma |A.4| the sum in ( |A.7| I 
is equal to 



E(-d j 



j=0 



j)\i + a-0-l 



= (-1)* 



a-N- 1 
/3-iV 



= (-!)'' 



(a-N - 1)! 
(a-^-l)!(/3-iV)!- 



By putting this into l |A.7| l, we see that the generic term of i writes 



/ 3!(i + a -iV)!(i + «-/?-!)! 
1 J + a)!(a-j8- 1)109 -JV)! 1 1 

which ends up the proof of Lemma |A.7[ 
Theorem A. 8. The matrix A~ 1 B is given by 



if 0\fi + a- -1\ fi + a 
NJ\ a-0-1 )'\ N 



□ 



A~ Y B = 



(-iy 



N f 0\ fa - + N - 1\(N - 1\ .fi + a 



i + a- 0\N 



N 



N 



0<i<N-l 
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Proof. Referring to Lemmas |A.5| and |A.7| we have that 

A~ l B = U(L~ 1 B) 



0<i,j<N-l 



'-"•(^)c:r/_r)/Cr 



0<i<N-l 



The generic term of A X B is 



The foregoing sum can be easily evaluated as follows: 



(A.8) 



N-l 

E 



A fj + a - p - 1 

a- (3-1 



i + a- 13-1 
a- 13- 1 



N-l 



JV-1 



i + a — [3 — 1 i 
a- (3-1 J ^ 





- a — 


P-l\ 




a — 


(3-1) 


:e 


+ a 






■f a 




/? + 


N- 




+ a 








j + a - /3 - 1 
i + a — /3 



TV 



j + a - (3 



a-(3 + N-l\(N-l 
N 



□ 



Putting this into I A.8 1 yields the matrix A 1 _B displayed in Theorem A.8 
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